Symmetry, Integrability and Geometry: Methods and Applications 



SIGIVIA 8 (2012), 009, 50 pages 



Lessons from Toy-Models for the Dynamics 
of Loop Quantum Gravity"^ 

Valentin BONZOM t and Alok LADDHA * 

^ Perimeter Institute for Theoretical Physics, 31 Caroline St. N, 
ON N2L 2Y5, Waterloo, Canada 
E-mail: vhonzom@perimeterinstitute. ca 

Institute for Gravitation and the Cosmos, Pennsylvania State University, 
University Park, PA 16802-6300, USA 
E-mail: alok@gravity.psu.edu 

Received October 11, 2011, in final form February 24, 2012; Published online March 07, 2012 
http://dx.doi.org/10.3842/SIGMA.2012.009 

Abstract. We review some approaches to the Hamiltonian dynamics of (loop) quantum 
gravity, the main issues being the regularization of the Hamiltonian and the continuum limit. 
First, Thiemann's definition of the quantum Hamiltonian is presented, and then more recent 
approaches. They are based on toy models which provide new insights into the difficulties 
and ambiguities faced in Thiemann's construction. The models we use are parametrized 
field theories, the topological BF model of which a special case is three-dimensional gravity 
which describes quantum flat space, and Regge lattice gravity. 
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1 Introduction 

The key feature of general relativity is the equivalence between gravity and geometry, so that 
geometry itself becomes a dynamical entity. Loop Quantum Gravity (LQG) [8, 87, 93, 105] 
is a quantization framework which particularly emphasizes those geometric content. It makes 
sense of quantum geometry of space, with quantum numbers of areas and volumes. And it is also 
more than that. It has been put to a rigorous mathematical status, which makes clear that it 
is a generic way to quantize background independent theories. This means that no background 
metric is available, and that general covariance must be satisfied at all steps. A key theorem 
asserts the uniqueness of the quantization map when diffeomorphism invariance is required [82]. 
As a result, the basic excitations are supported by embedded graphs, known as spin networks. 

After the completion of the kinematical framework, the task was to make sense of the Hamilto- 
nian of general relativity. This implies two steps: i) regularization within the loop quantization, 
and ii) the continuum limit. In a remarkable series of papers, [106, 107, 108], this task was 
carried out by Thiemann, and as a major success the quantization was shown to be anomaly- 
free, which means no classical degrees of freedom of the regularization are left over through 
quantization. 

This operator in addition to its technical precision has a number of interesting properties like 
diffeomorphism covariance. It has also led to detailed quantization of mini-superspace models 
like Loop Quantum Cosmology (LQC). Although its use in the current formulation of dynamics 
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in LQG is rather limited (spin foam models [10, 98] and alternative approaches like the Master 
constraint in the canonical theory [111] have been introduced to circumvent the difficulties faced 
by Thiemann's construction), as far as canonical LQG is concerned it still retains a fundamental 
importance. There are already a large number of reviews on this subject. Our aim in this review 
is not merely to repeat what has been said above, but instead to offer a frank account of the state 
of the art in the field, focusing on the most difficult aspects, and some of the different ways in 
which researchers have tried to go beyond Thiemann's construction. 

Quickly after the release of Thiemann's proposal, several criticisms were made, see for in- 
stance [64, 81, 104]. As usual in quantum theory, the construction is plagued by numerous 
ambiguities, but those papers also call into question the physical content of the continuum limit 
induced by the quantum Hamiltonian. In a nutshell, the difficulty with Thiemann's Hamilto- 
nian is that although the quantum algebra is mathematically anomaly- free, it was not possible 
to check explicitly that it implements the classical (Dirac) constraint algebra at the quantum 
level. In [104], this issue was traced back to the fact that the quantum Hamiltonian acts in an 
ultra-local way, only at the vicinity of nodes of the spin network excitations, so that physical 
states lack long range correlations. 

In spite of those criticisms, there has been very few modifications of the proposal (see [2] for 
a recent improved regularization) . However, as usual in physics, toy models are useful to get 
insights which should survive in the full theory, because they enable to disentangle the main 
issues. Here we will focus on 

• parametrized field theory in two dimensions, 

• topological BF field theory, 

specific case of the latter, 2+1 gravity, 

• a canonical approach to discrete gravity (Regge calculus). 

They aim at going deeper into the structure of the loop quantization, and explore new alterna- 
tives which are, sometimes surprisingly, fully amenable. Let us say a few words on those models 
(to avoid overflowing of references at this stage, we put them in the corresponding sections). 

Due to the infinite number of ambiguities of the standard quantum Hamiltonian, there is 
a vast room for improvement in the definition of the operator. In the context of parametrized 
field theory, we will argue that the off-shell closure of the constraint algebra is a key requirement 
which not only has the potential to reduce the number of ambiguities, but also that in simple 
models, this requirement leads to physically correct results. The main results appeared in 
[77, 78, 79]. 

Three-dimensional gravity is well-known as an exactly solvable topological field theory [113], 
where transition amplitudes have been evaluated [115]. However, one can try to forget as most 
as possible that it is topological, and treat it instead just as geometrodynamics with the loop 
quantization. The main question is then: What is the Wheeler-DeWitt equation, which would 
quantize the scalar Hamiltonian of the theory? Hence, this model addresses, in a successful 
way as we will show, the issue of the regularization in the loop quantization, and proves that it 
is possible to reproduce non-trivially the expected results of the topological formulation. The 
main references for the work we review are [27, 29, 115]. 

Three-dimensional gravity is a special case of the topological BF theory [25, 72] which can be 
defined in arbitrary dimensions. Though gravity and BF theory are not equivalent in dimension 
four and higher, since they differ by the presence and the lack of local degrees of freedom, it 
turns out that they are still related [41, 58, 60, 72]. One of these relations consists in formulating 
gravity as a BF model supplemented with constraints, and this is the main road to the spin foam 
quantization attempts to quantum gravity. If one wants to understand how spin foams may be 
able to circumvent the difficulties faced in the Hamiltonian quantization, one should therefore 
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start with a good understanding of the topological model. However, until recently, its loop 
(and spin foam) quantization was not so much developed. It is mainly known in the three- 
dimensional case, thanks to [59, 88], and progress in higher dimensions have only appeared very 
recently [34, 37]. A pedagogical review is [10] but it is now ten years old. We consider the 
section on the BF model of the present review as an update to this review, proving it is now 
much better understood. 

Since general relativity describes the dynamics of the geometry, it is interesting to see how 
the quantization of the 2+1 scalar constraint describes flat quantum geometry. This geometry 
is naturally discrete since it is supported by graphs, so that it makes contact with the canonical 
approach to discrete gravity, which looks for a Hamiltonian respecting the gauge symmetries, 
and hence the physical degrees of freedom of the full theory. Such a Hamiltonian can be reached 
through the so-called perfect discrete action, obtained after coarse-graining. The material of the 
review is borrowed from [13, 14]. 

The literature is quite rich and we have certainly left interesting works behind. For instance, 
we have not been able to cover most of the work done in certain mini-superspace models. 
Similarly, very little will be said on spin foam models, though parts of the work presented in 
the review were motivated by linking canonical LQG to spin foams, especially [29, 88] in 2+1 
dimensions. In higher dimensions, the most recent works are [1, 3, 33]. 

As we will show, the technology used in loop quantum gravity, especially in the case of 
quantum flat space, relies heavily on SU(2) re-coupling methods. While the basic tools are very 
well-known for a long time [112, 119], research in loop quantum gravity has led over the last years 
to interesting side-products on advanced quantum angular momentum theory. They include 
exact results, such as new recursion formulae for Wigner 6j-symbols and more complicated 
objects inspired by spin foam models [31, 33], and also new asymptotic expansion in the regime 
of large angular momenta. In particular, the expansion of the 6j-symbol is now understood 
beyond the usual leading order [32, 54, 55], asymptotics of arbitrary Wigner coefficients with 
some large and small quantum numbers have been obtained in [28] completing the exciting 
results of [116, 117, 118] based on an improved WKB approximation. The coherent states 
technology introduced in the context of quantum gravity in [83] has been very useful in deriving 
new asymptotics for 15j-symbols and other objects from spin foams [20, 21]. 

Such developments in quantum angular momentum theory have not appeared only in the 
quantum gravity community [5, 7], and is now considered as a key topic in modern physics [6], 
beyond the traditional use in spectroscopy and atomic/molecular physics, for instance in quan- 
tum computing [38], in the presence of topological order [80]. This is not so surprising since 
the reference model for topological order and fault-tolerant computing is the Kitaev model [75], 
which is actually a lattice BF model with the group Z2. We think those progress on quan- 
tum angular momenta in the study of quantum flat space is one of the major consequences of 
quantum gravity outside of the field and would deserve their own review. 

Such models lying at the border between toy models for condensed matter and for quantum 
gravity have been recently proposed in order to investigate the second key issue, the continuum 
limit in background independent quantization schemes. Since the project has been started very 
recently, we will not comment further, but refer to [16, 45] instead. 

A different road towards a generic framework for taking the continuum limit is under ac- 
tive development: tensor models and group field theories, which generalize matrix models to 
dimensions higher than two. A paragraph on those models is included but we refer to the 
reviews [57, 92] and more recently [69] for more details. 

Since we review quite advanced works, on the dynamics itself, we include for completeness 
some basic material of the loop quantization and its kinematical aspects in Section 2. Then, we 
give an introduction to Thiemann's construction in Section 3. In the interest of pedagogy, we 
restrict to the Euclidean Hamiltonian constraint. This is because the full Lorentzian constraint 
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is the sum of this EucUdean term with an additional piece, so that it is in any cases necessary 
to understand the Euclidean part. A small paragraph on the master constraint program which 
was launched by Thiemann after criticisms of its construction is included. Throughout this 
section, we deliberately focus on the unsatisfactory aspects of the construction, for two reasons: 
i) its mathematical consistency is already emphasized in previous reviews, ii) to stress the issues 
which we address later using toy models. 

While our presentation of the standard quantum Hamiltonian may sound a bit pessimistic, 
this is to contrast with the present situation which seems better to us. In the following sections, 
we move on to present some exciting results which have been obtained in the very last years. 
In Section 4, we consider a simple two-dimensional generally covariant field theory to extract 
some interesting lessons. In Section 5, the diffeomorphism constraint itself is shown to admit 
a quantization rather satisfactorily if one is ready to examine all possible choices we are afflicted 
with in LQG. In Section 6, we present a different approach based on the topological BF theory. 
It contains some review material on the BF model, together with the most recent quantization 
of the scalar Hamiltonian. A discussion on how to get beyond the topological model is also 
included. Some conclusions are presented in Section 7. 

2 Background material 

This section presents some background material on loop quantum gravity. More details can be 
found in the reviews [8, 93, 105]. 

Yang— Mills phase space. The (unreduced, i.e. before imposing the constraints which 
generate the dynamics) phase space of 3+1 general relativity in Ashtekar-Barbero variables 
is quite similar to the Yang-Mills phase space. On the canonical surface S, we have a SU(2) 
connection 1-form A^^, where i = 1, 2, 3 is a su(2) index and a is the form index, and the canonical 
momentum Ef, which satisfy the Poisson brackets 

The momentum is the non-Abelian version of the electric field. In 3+1 general relativity, it is 
the triad of the intrinsic geometry on S, related to the 3-metric qab by qq"''' = 6^^EfEj. Here 7 is 
a free parameter of the classical theory known as the Immirzi parameter. The Yang-Mills phase 
space is the same with 7 = 1. In three-dimensional Riemannian gravity, the phase space is also 
the same, with 7 = 1 (but the form index just takes two values when S is two-dimensional). 

Holonomies. From a connection, locally seen as 1-form on S with special rules for gauge 
transformations, one can build holonomies, whose geometric meaning is to enable parallel trans- 
port of vectors. The holonomy along a path £ : [0, 1] — )• S is Ui{A) = Vexp — J^A, where V 
denotes the path ordering. The loop quantization also features generalized connections, which 
are objects which assign holonomies to any paths. 

Kinematical states of the loop quantization. The loop quantization provides a natural 
kinematical Hilbert space to quantize the Yang-Mills phase space. The idea is to probe the 
connection only through a finite number of variables, similarly to a Wilson line which probes 
the connection along a given path. We introduce the holonomies Ue{A) of A along the path e. 
Given a graph F, with E edges and V vertices, and a function / over (SU(2))-^, we form the 
cylindrical function iprj S Cylp supported on F as 

i;rj{A) = f{UeAA),...,UeAA)). 

We are moreover interested in SU(2) gauge invariant states. Gauge transformations act on 
holonomies only on their endpoints. If /i is a map from S to SU(2), then the holonomy transforms 
as Ue{A^) = h{t{e))Ue{A)h{s{e))~^ , with t{e), s(e) being respectively the source and target 
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points of the path e. When focusing on a single graph F, this reduces gauge transformations 
to an action of SU(2)^ on the set of cylindrical functions over T. So from any function / over 
SU(2)^, one gets an invariant function by averaging over the SU(2)^ action. 

The algebra of such functions has a natural SU(2)^-invariant inner product which comes 
from the Haar measure on (SU(2))^, d/ir = He^ffe' giving 

{(pr,h\iJrj) = / Y[dgeh{gi,. . . ,gE)figi,. . . ,gE)- 

e=l 

By completing this space for the corresponding norm, one obtains the Hilbert space Tir = 

L2(SU(2)-^/SU(2)^,d^r). 

The union of all functions which are cylindrical with respect to some graph is denoted 
Cyl = Ur Cylp. It comes equipped with a natural Hermitian inner product inherited from 
the measures d^r- The completion of Cyl then leads to the kinematical Hilbert space of the 
loop quantization ^kin- The latter can be seen as the set of square-integrable functions over the 
set A of generalized connections, ^kin = L'^{A,dij.\i^), where /^al is known as the Ashtekar- 
Lewandowski measure, ^kin can also be constructed from the set of Hilbert spaces Tir for all 
possible graphs T, using projective techniques. Importantly, the measure //al is invariant under 
diffeomorphisms of S. 

The spin netw^ork basis. It is an orthogonal basis of the above L2(SU(2)^/ SU(2)^, d/^r) • 
The first idea is to expand any element of L^(SU(2)^, d/xr) onto the matrix elements of the 
irreducible representations of SU(2)^, labeled by E spins je G N/2. Then, rotation invariance at 
vertices imposes to contract the matrix elements of edges meeting at a vertex with intertwiners ty . 
If u is a node of F with ingoing links (ein), and outgoing links (eout), then an intertwiner is 
a map: ®ei^T~ije ~^ ®eout^>) (where Hj is the carrier space of the spin j representation) which 
commutes with the group action. Thus, to form the spin network basis on we need to 
assign spins to edges and a basis of intertwiners at each node. For a trivalent vertex with 
ingoing edges for example, there is a single intertwiner tj^jjia • ^ii ® '^02 ® "^js ^ C up 
to normalization, its components in the standard magnetic number basis being the Wigner 
3jm-symbol (ji, mi; j2, m2; js, msltjjjjja |0) = ms)- Hence, spin network states are the 

following SU(2)^-invariant functions 

E V 

Sr"''"Hgi,---,gE) = ^ Jj0e,?TT-e|5'e|ie,f^e) JJ('X'einie,"Ze|tt,| Igieout Je,?^e>, 
{me,ne}e=l v=l 

which form an orthogonal set. This enables a transform between gauge invariant functions on F 
and functions over the colorings {je,Ly). 

Phase space on a graph. There is a natural classical phase space on F, inherited from 
the full phase space, whose quantization leads to L^(SU(2)'^/ SU(2)^, d/ir)- Each link e carries 
a group element ge S SU(2), which can be seen as the holonomy of the connection. In this 
picture, each node carries a local frame, and ge performs parallel transport from the frame of 
the source node of e to its target node. The momenta Ef give rise to fluxes when integrated 
along surfaces of codimension 1 on S. If it is done infinitesimally close to a node, source of 
a link e, we get the flux variable El G su(2), defined with respect to the frame over the source 
node. The Poisson bracket is the structure of the cotangent bundle over SU(2)^, 

{El,Ei}=e'\Et {Ei,ge}=geT\ 

where the matrices (Tj)j=i.2,3 are anti-hermitian matrices^ satisfying [r*,r-^] = e\T^ . All other 
brackets vanish. One can obviously define a flux variable in the frame of the target node of e, 
using the adjoint action of the holonomy on the initial flux, E^. = Ad((7e) E^- 

^They read: = — f in terms of the Pauli matrices ((t'). 
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Figure 1. The dual graph, in dashed Unes, which supports the phase space, and the triangulation. The 
fluxes are the normals to the edges of the triangle, they encode the local embedding into flat 3-space. 
Their dot products are proportional to the cosines of the angles ipee' ■ 



Geometric content. For simplicity, we consider T to be dual to a triangulation of S. Fluxes, 
which come from the triad, encode the intrinsic geometry of the triangulation. The main reason 
is that the Gaufi law, which generates SU(2) transformations at a node v, reads 

^ El- ^ El = 0. (2.1) 

e outgoing e ingoing 

Thinking of the node as dual to a tetrahedron, and the links dual its triangles, each of them 
carries a flux and this allows to interpret the above relation as the closure of the tetrahedron, 
where fluxes are the (oriented) normals to the triangles (in the frame of the tetrahedron) . Areas 
are given by their norms, and dihedral angles by their dot products. This way, the geometry of 
each tetrahedron is completely specified. However, different tetrahedra will generically assign 
different lengths to their common edges. Such geometries are described in [47, 48, 61, 100] and 
are known as twisted geometries. 

A similar interpretation holds in 2+1 gravity. Since fluxes are still three-dimensional vectors, 
one should rather think of nodes as dual to triangles locally embedded in flat 3-space. The 
GauB law on a node is the closure of the dual triangle, so that fluxes based on a node represent 
the normals to the dual edges in the frame of the dual triangle, as shown in Fig. 1. Norms and 
dot products of fluxes give access to the lengths and 2d angles within the triangles. The closure 
relation in particular asserts that 

COS0e.e.= "|t. " (2.2) 

which is the standard formula to compute the angles of flat triangles from their lengths, namely 
the law of cosines. 

In contrast with the 3+1 case, the geometries encoded into fluxes in 2+1 dimensions assign 
lengths unambiguously to the whole triangulation and are thus Regge geometries. 

It is also possible to introduce dihedral angles between the top-simplices of the triangulation, 
which give a notion of extrinsic geometry, using fluxes and holonomies. The Hamiltonian of 
general relativity is precisely a relation between extrinsic and intrinsic geometries, and will be 
described later in the flat case. 



3 Quantum Hamiltonian constraint: Thiemann's definition 

In this section we review the quantization of Hamiltonian constraint, as proposed by Thiemann 
in [106, 107]. Instead of focusing on the details, we focus on the basic underlying ideas and what 
in our opinion are its major shortcomings. 
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In a remarkable series of papers, Thiemann proposed to quantize the Lorentzian Hamiltonian 
constraint of canonical gravity by using in a crucial way certain relations involving Poisson brack- 
ets at the classical level, together with suitable choices of intermediate regularization schemes. 
The final Hamiltonian constraint turned out to be well defined on ^kin- It showed glimpses 
of generating an anomaly free quantum Dirac algebra, as the commutator of two Hamiltonian 
constraints, when acting on any spin network state Tg in Hkim gave a null state when averaged 
over spatial diffeomorphisms. 

However as we will see below, the definition given by Thiemann is unsatisfactory for several 
reasons. This is not a criticism of the construction itself as it was the first ever rigorous definition 
of the Hamiltonian constraint, but as it has been over fifteen years since Thiemann's seminal 
work, it is perhaps good to take stock of the current state of the art in the field. 

As we argue below, the most worrisome aspect of Hamiltonian constraint is that there are 
infinitely many ambiguities which feed into its definition. Each such choice could give rise to (in 
the continuum limit) an incquivalcnt quantization of the Hamiltonian constraint. The preferred 
choices currently used in the literature are based on subjective criteria rather then appealing to 
certain physical requirements like local Lorentz invariance. 

The classical, density one, Euclidean Hamiltonian constraint is given by 

r F''' E°-E^ 

He[N] = / d\N{x)et i- 

(we have set G = 1). As we will see shortly, it is only for density weight one, that the continuum 
constraint operator can be well defined on 7^ kin- However it is important to note that, it is 
a priori not required that the domain for the quantum Hamiltonian should be in T^kin, and 
could very well be a subspace of Cyl* in which case other density weighted constraints could 
also admit a continuum limit. 

One of the ingenious features of Thiemann's construction is the following. The density one 
constraint, although non-polynomial in the basic fields could be written as 

He[N] = [ N{x) Tt{F{x) a {A{x), V{R)}), 
Jy, 

where R is any open region containing the point x. V{R) is the volume of this region, as 
a function of the field E. 

As the elementary variables in the theory are holonomies along (piecewise analytic) edges 
and fluxes along (semi-analytic) surfaces, we need to first approximate He[N] by a sequence 
of functionals of holonomies and fluxes such that the limit of the sequence (in the suitable 
sense) is He[N]. This is accomplished as follows. Choose a cell decomposition of the under- 
lying manifold S (it can be simplicial or cubical for convenience). For each 3-cell A in T, 
choose a vertex v(A), a family of edges Se = {si(A), . . . , sMs('^)} and a family of loops 
Sl = {ai(A), . . . , Q:Mi(A)} and consider 

Ms Ml 
AeT(e) i=l J=l 

where C*"^ are constants independent of e. Each choice of (T'(e), {C^'^},Se,SL ) such that, 

lim HT^,)E[m = HE[N] 

is an acceptable net of regulators. Classically one choice is as good as the next (in the limit of 
infinite refinement) but quantum mechanically it is not clear which of the choices are indistin- 
guishable in the continuum limit. 
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As i?j'(e)[A^] is a function of holonomies and volume, it can be quantized on ?^kin- Given 
a spin network s = (7,j,c) let T(e,7) be a one parameter family of triangulations adapted to 
the graph 7. By adapted all we mean is that each vertex v' in V{'y) should be either in the 
interior of precisely one cell A G T{e, 7) V e, or v' be a vertex of T(e, 7) such that for all e > 0, the 
number of 3-cells saturating v' remain constant. One can now prescribe the action of Hx(e,'y)[^] 
on |s) as, 

Ms Ml 

Hne,,m\^)= E EE^'''Tr(/i„,(^)/.,^(^)[y^),y(i?,(^))])is). (s.i) 

AGT(e,7) i=l J=l 

Let Ti^ be the subspace of Tikm spanned by spin networks based on graph 7. There are two 
definitions of volume operators available in the literature. One is known as AL (Ashtekar- 
Lewandowski) volume operator and the other is known as RS (Rovelli-Smolin) volume operator. 
We will restrict ourselves to the AL volume operator in this review^. The equation (3.1) describes 
the action of iJ-p^^ [N] on H-y for any 7. Whence we have a densely defined operator Hrpf^^^^-^ [N] 
on the (dense domain of) entire Hilbert space 

So far the parameter space regularizing the Hamiltonian constraint is uncountably infinite 
dimensional. This is by any account quite huge. This itself is not an issue because it might 
happen that different regularization schemes lead to same continuum operator (or at least the 
the same physical spectrum of the theory) . 

One would also expect that putting physical requirements like diffeomorphism covariance 
will restrict the choice of triangulation. It turns out that the requirement of diffeomorphism 
covariance places the following restrictions on the allowed choices of triangulations. Given a pair 
of quintuples (7, f , A G T^, 5e, Sl) and (7', v' , A' £ T[e'),S'^,S'i) such that (7, v) is diffeomorphic 
to (7', f'), we require the entire quintuples to be related by some diffeomorphism for any choice 
of e, e'. 

Let us introduce U{(j)) as the unitary operator corresponding to a diffeomorphism (j) which acts 
on spin network states like U{(j))Ts = T^-i^. It can be easily shown that the above requirements 
imply that V0 G Diff(S) 

U{ct>)HT(^^^,)[N]U{^-^) = Ui<p')HT(^-.^-y),e)[rN], (3.2) 

where </)' is a diffeomorphism not necessarily equal to (p. As we will see in the next section, this 
would imply that the continuum Hamiltonian constraint in diffeomorphism covariant. 

However the covariance requirement, though non-trivial, does not reduce the parameter space 
of allowed regulators significantly^. 

We emphasize that this is the only physical requirement one places on regularization. 

3.1 Continuum limit 

There are various operator topologies in which one can examine the convergence (or lack thereof) 
of Hj'(^^^y-)[N]. The finest topology used in the literature so far is the so called URST topology^, 
in which one defines the continuum limit on ?^kin itself. This is quite an important construction 
so let us review it a bit more. 

^If one insists on defining continuum limit of Ht(S) [N\ on 'Hkin then it is essential to use AL volume operator. 
However if one instead defines continuum limit on some subspace of Cyl* then as shown in [81], RS volume 
operator is a perfectly viable choice. 

^Although it should be emphasized that to the best of our knowledge the only explicit construction of a co- 
variant family of regulators available in the literature is given in [106], but that of course does not mean that it 
is the only choice, in fact as explained in [8], there are infinitely many. 

''URST is an abbreviation for Uniform by Rovelli, Smolin and Thiemann. 
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One corollary of requiring the one parameter family of regulators Re, 7 to be covariant is that 
for any two e, e' > and for any vertex v G ^^(7) the corresponding family of loops 5^,^,(6), 
Sl,v{^') are diffeomorphism covariant. This implies that 

(^|^r(,)[iV]|s) = (vI/|i7^(,,)[iV]|s) 

G ^diff- One can now define the continuum limit of the Hamiltonian constraint (on ^kin) 
the following way. H[N] is a continuum limit of Hj'i^^-^\N] if, given any (5 > 0, there exists 
eo((5) > such that for all e < eo 

|(*|i^[iV]|s) - (vI.|^r(,,^(,))[iV]|s)| < 5. (3.3) 

This is trivially true because any H(^^^-^[N]\s) and H f^^^-^lNWs) are related by a diffeomorphism, 
so that (^'|lfy(e,7(s))[-^]|s) is actually independent of e. 

Whence in the above mentioned topology, Hrp^^^-^^^N] is already a continuum operator for 
any eo. Two properties of this construction are worth noting. 

1. Since eo is independent of (|s), ^ and N), the convergence is uniform. 

2. If two operators when acting on any state in ^kin give diffeomorphically related states, 
they are identical. This in conjunction with the equation (3.2) implies that the continuum 
operator is diffeomorphism covariant! 

If we do not care about demanding uniform convergence (and it is not clear to us if there is any 
physical reason to care), then we can define the continuum limit on infinitely many subspaces 
of Cyl*. That is, if we define continuum limit via, 

[^[iV]'M/]|s) = lim^[^r(e)[A^]|s)], 

where the convergence is only point wise (the notion of continuity used is same as the one in (3.3) 
except that eo = eo((5, ^',s)) then on a large class of states in Cyl* H[N]' could be well defined. 
Let us illustrate this point using a following schematic example. Let 

where A is some set of spin networks {A can be uncountable) and / are functions of the spin 
networks defined on a suitable domain, (e.g. / could be functions of vertices of s', or they could 
be functions of vertices and germs of edges incident on those vertices, or they could be functions 
of contractible loops contained in s' etc.). We refer to / as vertex functions^. Schematically 

m 

i)ev(7(s)) 1=1 

where due to diffeomorphism covariance of regulators, the coefficients a^j are independent of e 
and |s^j) = i/((/)e)|sJ5j) for some reference spin network |sj]j). 

Now let the set A be such that if it contains some state then it necessarily contains all its 
diffeomorphic images. Furthermore let us for the sake of illustration assume that there exists 
precisely one vertex vq for which |s^^j-) G A. Whence 

m 

t)ev(7(s)) 1=1 



^See [64, 81] for more details. 
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m 



m 



v&V{-y{s)) 1=1 




Thus 



m 



i?[Ar]'^r|s) = N{vo)Y 



1=1 



Such subspaces of Cyl* on which continuum hmit of the Hamiltonian constraint is well defined 
are known as the habitats. Our schematic computation illustrates that potentially there are 
infinitely many such habitats. One explicit construction of a habitat is given in the seminal 
paper [81]. Essentially it is characterized by the property that the function / above is a (smooth) 
functions of vertices of the s. This habitat turns out to be a sufficiently small extension of 'HdHi- 
We will have many occasions to use this habitat in this review and henceforth refer to it as 
Lewandowski-Marolf (LM) habitat. 

This finishes the summary of the construction of candidate Hamiltonian constraints in LQG. 
Let us, for the benefit of the reader summarize some of the worrying aspects of the construction. 

1. The parameter space underlying allowed choices of regulators is infinite dimensional. This 
would not be a problem if the continuum limit of the Hamiltonian constraint is independent of 
such choices, however this is far from being the case, e.g. in the URST topology, each pair of 
choices which is not diffeomorphically related to each other will give rise to a different operator. 

2. It is a priori not clear what the domain of the Hamiltonian should be. It could either be 
a dense subspace of kinematical Hilbert space, or a habitat. 

3. There are further ambiguities. In all the constructions available in the literature so far. 
Lapse is a multiplicative factor. However one could conceive type of regularizations where the 
lapse function is included in the definition of the plaquette one uses to define curvature [79]. 

4. There is a huge amount of choice available in defining curvature operators. One such 
choice in particular has received some attention in the literature. This has to do with the 
representation underlying the curvature approximant in the Hamiltonian constraint [93, 94]. 
The simplest choice for such a curvature approximant is 



where Tij is the trace in some representation of spin j, and dj is the appropriate normalization 
factor. In fact as in LQG one often employs "state dependent" regularization schemes, this 
ambiguity could be generalized further where e.g. 



Here Trj^[he]{A) is the gauge- variant spin network state based on a single edge e, and v £ e. 
whence the choice of representations used to define curvature approximants could themselves 
be state dependent (and in this sense dynamical). 

5. Yet another ambiguity which seemed to have gone un- noticed in the literature until the 
arrival of improved dynamics in LQC is the following. Notice that the flux operator when acting 
on a spin network only depends on the intersections of the edges of the spin networks with the 




where the trace is in the fundamental (j = 2) representation. One could instead choose approx- 
imant to be 
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underlying surface. That is, even if we scale the coordinate area of the surface arbitrarily, the 
action of flux operator on a spin network does not change as long as the intersection numbers 
remains the same. This means the following. Consider latticized classical curvature is given by 

However one could equally well consider a different approximant of the type 

where Aab is a plaquette based at v of co-ordinate area e^, and is a surface also of co-ordinate 
area e^. Classically the second term is higher order in e than the leading order term, however 
quantum mechanically the difference between states obtained by action of two approximants on 
a given spin network state will be a finite norm state. Whence this is a genuine ambiguity in 
the choice of curvature which can drastically affect the definition of the Hamiltonian constraint. 
In essence, it is this choice which was exploited in [9] to define a hamiltonian constraint in LQC 
which gave rise to better semi-classical properties of LQC. 

Thus the situation may look rather worrying. We have infinitely many possible definitions 
of continuum Hamiltonian constraint and it is a priori not clear which if any of these choices are 
physically preferred (better semi-classical properties). 

It should be mentioned that as explained in [105], it is conceivable that as the final spectrum 
of the theory consists of diffeomorphism invariant states, at least some of the choices which 
give rise to the kernel of the Hamiltonian constraint that are diffeomorphically related may not 
matter. Although this argument is intuitive and to the best of our knowledge not analyzed in 
detail, it is an important counter-point to the grim view that we have presented above. 

Finally a worrisome feature of H[N] not necessarily related to the infinitude of ambiguities 
is its ultra-local character (see [93, 104] for details). Consider a spin network state |s) and let 
have support in the neighborhood of a vertex vq G V{'~f{s)). Then H[N] (assuming it is defined 
on 'Hkin) will give rise to a linear combination of states each of which is "different" from s only 
in the neighborhood of vq (which means only the edges, which are incident on vq will undergo 
any change under H[N] action.) This statement would remain true even if we worked with 
arbitrarily fine graphs and hence even when we consider say (kinematical) semi-classical states 
based on a single graph. Naively we expect the Quantum constraint equation in this case to 
"approximate" classical constraint equation which are elliptic and hence not ultra-local. Thus 
the semi-classical limit of such Hamiltonian constraint might turn out to be incorrect. One 
way around this issue suggested in [104] is to tailor the choice of loops underlying the curvature 
approximants in such a way that it intersects many vertices in a given graph. Such a construction 
has been done independently in the context of BF theory for 2+1 gravity, detailed in Section 6. 
Finally the criticism regarding ultra-locality has itself been criticized in [105] on the grounds 
that non-locality might emerge when one considers physical states which live in Cyl*. 

3.2 Action of H[N] on ^kin 

In the URST topology we have a continuum operator on T^kinj 

iV(s) 

H[N]\s) = HT^,^^[N]\s) = ^aj\s,), 

1=1 

where sq is fixed once and for all^. Each s/ has two more vertices and one additional edge than s. 
The two new vertices that are created are planar, and are in the kernel of AL Volume operator. 

^There is a minor subtlety here. As the HT(to)[-^] when acting on |s) must be finer then the graph 7(s). 
eo implicitly depends on 7(s). 
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This feature plays a key role in the "anomaly free" nature of the Hamiltonian constraint, when 
it is defined on 'Hkin- 

3.3 The commutator and its unexpected triviality 

One of the interesting properties of the Quantum Hamiltonian constraint defined by Thiemann 
was its "on-shell closure". That is the commutator of two continuum Hamiltonian constraints 
vanished. Although consistent, a closer look at the computation shows that the underlying 
mechanism for this commutativity is not too restrictive (as far as the regularization choices are 
concerned) and does not rely on any deep or non-trivial structural aspects of the Hamiltonian 
constraint. 

This is because the commutativity only relies on two features of the regularization. 

1. The density one nature of the constraints ensure that the corresponding operators are 
"ultra- local" (whether we define them on T^kin or on Cyl*) and can not give rise to one 
forms A^VqM — MVaN in the commutators. 

2. The diffeomorphism covariance condition placed on the regularization scheme (which is 
the only physical requirement fed into the quantization so far) guarantees that the commu- 
tators at finite triangulation is anti-symmetric in the lapses. This can be seen as follows. 
As before, we will merely sketch the proof. Details can be found in [81, 105]. 

Let us schematically write the Hamiltonian constraint operator at finite triangulation as 

^T(e)[iV]|s)= N{v)d,{v)\s). 

veV{-y) 

Here O^iv) is a non-local operator which depends on a set of loops (based at v), set of edges 
and the Volume operator. 

The essence of diffeomorphism covariance of the regularization scheme is the statement that 
for any graph 7 and e, e', 3 a diffeomorphism (jje^t'il) which preserves 7 and 

a(t')|s) =?7(0,,,,,^)4'(^)|s) (3.4) 

Vv G ^(7). Consider two lapse functions A^, M and as the Hamiltonian constraint operator 
is linear in lapse we can (without loss of generality) assume that given a graph 'j N, M have 
support only in the neighborhood of vertices uat, vm respectively. Clearly the commutator 
between Hj'(^^-^[N], Hrp^^,-^[M] is zero for all e, e' if vn = vmJ whence the continuum limit of 
such a commutator (in URST topology on ^kin or some weaker topology on a suitable habitat) 
will be zero. 

If vn 7^ Vm then it can be easily shown that (essentially using the covariance of regulator as 
in (3.4)), 

[HTie)[N]HTie')[M] - N ^ M] |s) = N{vN)M{vM){U{H^,e')) - l)d,{vN)d,>{vM)\s). 
Clearly in the URST, this operator vanishes as 

mu{cp) -i)\s) = o 

V(^'| G "Hdiff, (f) € Diff(S), |s) G "Hkin- The argument can be extended for the continuum 
Hamiltonian constraint defined on Lewandowski-Marolf habitat (see [81] for details). Although 

'^This statement is true if we use AL volume operator. However as shown in [81], on the LM habitat, even if 
we use RS volume operator, even though this statement will no longer hold, in the continuum limit [//[A'^], i^fAf]] 
will vanish if wjv — vm- 
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plausible, we are not aware of any arguments showing that this result will continue to hold when 
the domain of the continuum Hamiltonian constraint is any suitable subspace of Cyl*. 

If one works on the LM habitat, this commutativity is puzzling at first. This is because 
states in the LM habitat have the following properties. 

1. It carried a representation of the algebra of spatial diffeomorphisms. 

2. By its very construction it contained the diffeomorphism invariant states. 

Whence not all the states in the LM habitat are diffeomorphism invariant, and we do not expect 
the commutator of two Hamiltonian constraints to vanish on the habitat. In [64] it was shown 
that this apparent inconsistency could be avoided by a specific quantization of r.h.s. (Essentially 
what was shown in [64] was that there exists a quantization of q"'^ which gave an identically zero 
operator on the habitat.) This happens because of the following simple reasons. Consider the 
r.h.s. of the Dirac algebra (for density one Hamiltonian constraint) 

r.h.s. = / d^xq"'' {NVaM - MVaN) F(diff);,(x). 
Js 

It turns out that one can quantize q'^^Hb using identities introduced by Thiemann and the final 
operator at finite triangulation has the following form, 

T.hl^^,) = Y^HNVaM - MV.iV) ,d^(,)(t;(A)), 
AeT 

where 0^(e)(w(A)) is a composite operator built out of holonomies and fluxes and hence it is 
well defined on ^kin- Thus if the limit e — )• exists on some space (whether on 'Hkin or on LM 
habitat or some other subspace of Cyl*) then it will trivially be zero. 

Naively one would think that the chief culprit in the game is density one character of the 
constraint. That is, if we were to work with higher density constraints then at finite triangulation 
the operators will be of the type 

HTi4N]\s) = ^ E N{v)d,{v)\s), 
i,ey(7{s)) 

where we have placed tilde on H to emphasize on higher density nature, o is a positive integer 
which depends on the density of H. 

Due to the presence of e in the denominator, it is certainly conceivable that commutator 
of two such Hamiltonians could give rise to factors of the type ^(")^^(^+'^)~^^(^)^(^+'^) _ jf 
limit e — 7- exists on some habitat, then the continuum limit of the commutator would contain 
factors of the type MVaN — NVaM. Whence higher density constraints have a chance to yield 
a faithful closure. However as argued in [64], this is not enough. Unless the new vertices created 
by the action of Hamiltonian are such that, they can themselves admit a non-trivial action 

of -ff[A^], the commutator would never produce infinitesimal diffeomorphisms. 

Thus the necessary conditions to obtain off-shell closure are, higher density weight of the 
classical constraint, and a choice of regularization which produces "non-trivial" vertices in the 
sense explained above. So far, no quantum Hamiltonian constraint satisfying these conditions 
has been found. 

One could question the importance of off-shell closure in the formulation of LQG. After all, 
we have at our disposal the diffeomorphism invariant Hilbert space, and working within this 
arena, off-shell closure condition is super ceded. However the reasons for requiring off-shell 
closure are two fold. 
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1. As we argued above, the allowed choices of regulators which yield well defined continuum 
Hamiltonian constraint on either 'Hkin or some habitat are infinite in number. It is im- 
portant to reduce this ambiguity, not by appealing to subjective criteria, but to certain 
physical conditions. A key physical requirement would be local Lorentz invariance, which 
in canonical gravity is encoded in the Dirac algebra. Prom this point of view, demanding 
Dirac algebra to be faithfully represented in the quantum theory seems very appealing. 

2. Moreover, the fact that commutator continues to vanish even on non diffeomorphism 
invariant states seems like an extremely unpleasant feature of the operator. The fact that 
one can choose to quantize r.h.s. in such a way that it also annihilates such states does 
not seem very comforting as it is tantamount to the inverse metric being zero operator. 

It is sometimes argued, [105], that the LM-habitat is too small an extension of I^^jff to see any 
non-triviality or anomaly (I'^iff is a dense subspace of ^^iff which contains the diffeomorphism 
invariant distributions). However this is not true, as has been explicitly demonstrated in simpler 
contexts like Parametrized Pield Theories (PPT) [77]. We will summarize results obtained in 
this simple model below in Section 4. 

Whence as far as quantization of the Hamiltonian constraint in LQG is concerned, this is 
where things stand. One has a slue of regularization ambiguities each of which can potentially 
give rise to a different continuum operator. Non-trivial checks like the off-shell closure of the 
constraint algebra do not seem to go through and in simple models like Loop quantum cosmology 
and PPT, it has been shown that the regularization choices used in the full theory so far leads 
to either incorrect (in PPT case) or physically unsatisfactory (in LQC case) results. 

One of the key lessons that LQC teaches us with regards to LQG is the need to write the 
curvature operator as a concomitant of holonomy and flux operators. Even though the reason 
as to why this curvature operator is selected in LQC comes from detailed semi-classical analysis 
of physical observables, it is interesting to note that the curvature operator used in that case 
does not "fall into" class of regularizations considered for curvature in the full theory. 

We can now pose a following question. Is there a subset of choice of regularization ambiguities 
(choice of triangulation, choice of curvature, choice of operator ordering etc.) and a habitat 
in Cyl* such that quantization of (certain density weighted) Hamiltonian constraint and certain 
quantization of the diffeomorphism constraint gives a faithful representation of the Dirac algebra 
on the habitat. If we managed to find such regularization scheme, we could further ask if the 
choice of curvature operator had any (at least conceptual) similarity with the curvature operator 
used in LQC. 

3.4 The master constraint program 

Before going into the toy models at the core of the review, let us say a few words on the master 
constraint program. 

Criticisms on Thiemann's construction appeared quite quickly after publication of his work. 
But due to the difficulty of the challenge, few modifications aiming at improving it have been 
proposed, and efforts have instead concentrated on developing the path integral spin foam for- 
malism. In [111], Thiemann himself introduced the master constraint program, to revive the 
interest into the canonical quantization. He traced back the difficulties in the well-known pe- 
culiarities of the Dirac constraint algebra. In the latter, the scalar constraint does not generate 
a subalgebra, is not spatially diffeomorphism invariant, and the whole algebra closes with struc- 
ture functions which are not polynomials in the metric components. 

Hence, Thiemann introduced in [111] a re- writing of the scalar constraint as 
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The density-weight one Hamiltonian H{x) is squared, rescaled to y^det{q{x)) to make it back 
density- wight one, and finally integrated over the canonical surface S. That stands for the 
infinite number of constraints H{x) = 0. M commutes with the generator of spatial diffeomor- 
phisms. 

The price to pay for getting a single master equation is that a (Dirac) observable O is now 
defined by a non-linear equation, 

{0,{0,M}}|M=o = 0. (3.5) 

The program has been successfully tested in various situations, [49, 50, 51, 52, 53] and an 
explicit comparison between the Dirac and path integral quantizations is presented in [71]. The 
loop quantum gravity case is developed in [110]. 

At the technical level, the master constraint program suffers from ambiguities just like the 
previous attempt. At the conceptual level, its point is to avoid the requirement of reproducing 
the Dirac algebra. It is appealing since the diffeomorphism algebra indeed encodes gauge, 
non-physical degrees of freedom. However, the generators of the algebra also form the full 
Hamiltonian of general relativity, which means that solving the dynamics is equivalent to finding 
the equivalence classes determined by different gauge orbits. It seems to us that this intuition 
becomes more complicated in the master constraint program, in particular since the equation for 
Dirac observables is not anymore invariance along the gauge orbits but (3.5) instead. Obviously, 
because the master constraint program is a re-writing, the usual notions from gauge systems 
have some equivalent in this formulation, but it is not clear to us whether hiding the difficulty 
of working with the Dirac algebra is fully satisfactory. 

4 Parametrized field theory as an ideal toy model 

The main question coming from our review of Thiemann's construction is whether there are some 
choices of regularization and a habitat such that the Dirac algebra can be faithfully represented. 
Needless to say, posing this question in LQG is going to be pretty hard, so as a warm-up we can 
pose it in some simplified model. The model has to be be a field theory so that it has spatial 
diffeomorphism covariance, and it has to be generally covariant, so that in the canonical theory 
one has Dirac algebra as it's constraint algebra. 

Parametrized Field Theories (PFT) present us with just such a class of models. These are 
just free field theories on flat spacetime in d > 2 dimensions written in generally covariant guise 
by adding "fake" (pure gauge) degrees of freedom. For details we refer the reader to [70, 76]. 

Two dimensional parametrized massless field theories on a cylinder are the simplest toy 
models for canonical gravity as the constraint algebra just becomes the well known Witt algebra 
(Diff(5^) ® Diff(S'^)), and hence is a true Lie algebra. The corresponding gauge group (which 
would not even be defined as an ordinary Lie group [26]) in canonical gravity is ^pft = Diff(5'^) x 
Diff(S'^). This is very nice as in LQG we know how to handle spatial diffeomorphisms very well. 
Thus once we quantize the kinematical Phase space of the theory using background independent 
techniques that LQG provides for us, we can average over the entire gauge group of the theory 
^PFT using the Refined Algebraic Quantization method as outlined in [66]. So on one hand we 
have Hamiltonian constraint and diffeomorphism constraint just as in canonical gravity, and on 
the other hand we could complete the entire quantization program without worrying about the 
constraints and obtain physical Hilbert space of the theory. So we could pose the problem we 
posed above for LQG Hamiltonian constraint with a further consistency check. If we instead of 
averaging over ^pfTj quantized Hamiltonian and diffeomorphism constraints, is there a choice 
of their quantization such that the kernel of the constraints yield the physical Hilbert space 
obtained via group averaging techniques. With this motivating monologue out of the way, we 
delve into more details. 
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4.1 Brief review of polymer PFT 

Consider the following canonical field theory on a Lorentzian Cylinder {A4 = x R, 77). Cauchy 
slices are oriented circles coordinatized by the angular coordinate x S [0, 27r], with the direction 
of angular increase agreeing with the orientation of the circle. The theory has the following sets 
of fields. Two scalar fields , X~ which are known as embedding fields and their conjugate 
momenta, and a matter scalar field / and it's conjugate momenta^. 
The Structure of the model can be summarized as follows: 

Canonically conjugate embedding variables: {X^ {x),Il-^-{x)), (X~ (x), n_ (x)), X^{2Tr) 
= X±(0) ± 2tt. 

Matter variables: y='=(x) := -Kfif, with the brackets {f{x),TTf{y)} = 6{x,y), {Y^,Y^} = 
0, {Y^{x),Y^{y)} = ±{d.S{x,y)~dy5{y,x)). 

Density weight 2 constraints: H±{x) = [n-|-(x)X='= (x) it ^5^='=(x)^]. 

The constraint algebra is isomorphic to the Lie Algebra of vector fields on the circle. This 
algebra is more commonly known as the Witt algebra. 

Instead of working with the above constraints which form a true Lie algebra, we can also 
work with the "traditional" diffeomorphism and Hamiltonian constraints which form the Dirac 
algebra with structure functions. 

Diffeomorphism constraint: Cdiff generates spatial diffeomorphisms 

Cdiff (x) = H+ + H^ = [n+(x)X+'(x) + n_(x)X-'(x) + ^/(x)/'(x)] . 
Hamiltonian constraint: Cham generates evolution normal to the Cauchy slice, 

= ^ =\U+{x)X+'{x)-U^{x)X-'{x) + \(TTj + f'^)]. (4.1) 

y/X+'{x)X-'{x) ^ ^ ^ ^ ^ ' V ; 4V / •/ ;j V ; 

Constraint algebra: The Poisson algebra generated by Cdiff and Cham is the Dirac algebra, 

{Cdiff [TV], Cdiff = Cdiff [iV, M], {Cdiff [iV], Cham [M]} = Cham[i^AfM], 

{Cham[iV], Cham[M]} = Cdiff [/3(Ar, Af)], (4.2) 

wherein A^, M are shift vectors, A^, M are lapse functions and the structure function f3'^{N, M) := 
q''^{NVhM - MVbN) in (4.2) is defined by the induced spatial metric qab, 

qabdx^dx^ = -X+'X-'{dxf. 

Whence working with the density one constraints, we are led to a two dimensional avatar or 
constraint algebra of canonical gravity. On the other hand, working with density tow constraints, 
one obtains an algebra without the structure functions. The Lie algebraic nature of algebra 
generated by density two constraints is special to two dimensions. This happens essentially 
because in two dimensions (with one dimensional Cauchy slices) , the structure function changes 
from q""^ (for density one constraints) to qq""^ for density two constraints, but qq"''^ = 1. What 
we would like to check is if either of these two (or in principle any density constraints) admit an 
off-shell closure in quantum theory. 



*We refrain from introducing the Lagrangian formulation of the theory from which the canonical structure is 
derived, and from where the fields receive their titles (embedding fields and matter fields), as it is not relevant 
for our purpose. It suffices to note that we have a field theory with constraint algebra which is isomorphic to 
Dirac algebra in two dimensions and which has no true Hamiltonian. 
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4.2 Quantum theory 

A charge network s is a finite collection, 7(s), of colored, non-overlapping (except at vertices) 
edges, e, which span the range of the angular coordinate x (i.e. [0, 2it]), the colors being referred 
to as charges, and the collection of edges being referred to as a graph. Charge network states 
are in correspondence with charge networks and constitute an orthonormal basis similar to spin 
network states in LQG. 

4.2.1 Embedding sector 

• Charge network: = {j{s^), {k\, . . . ,k\)} where k\ are embedding charges whose 

range is specified by k\ G itk'^i ^ ~ l)---;'^^- Here A is a fixed, positive, 

integer-valued Barbero-Immirizi-like parameter. It is useful to define the 'minimum length 



increment', as a : = 



2nL 
A ■ 



Elementary variables: X^{x), Tg± [Il±] := exp [ — i Yl /e± n±] • 

e67(s+) 

Representation: T^± denotes an embedding charge network state. X^{x),Tg± denote the 
operators corresponding to the classical quantities X^{x), Tg±[Il±]. Their action is given 
by 



T,±T + 



and 



X^{x)Ts± := A^_^±r^±, 



where, for 7(5^) with edges 



h k\ 



if X G Interior (e^±) for 1 < < n^. 



\ h (k\ +k\ ) if X e et n efr^^,, for 1 < /± < (n± - 1), 



i/,(fc± ^^^ + k%) ifx = 0, 

if X = 2tt. 



'1 



\n[k^,±^t + k% 



± 



4.2.2 Matter sector 



Charge-network: s = {7(5 ),(/ ±)}, with the technical condition Yl ^ ± ~ ^ 

and G eZ for all / = 1, . . . , n^. Here e is a fixed real positive parameter with dimensions 
{ML)^2, e is also a Barbero-Immirizi-like parameter. 
Elementary variables: W^±^ [Y^] = exp[i X]e±e7(s±) C± /e± ^"^l- 

Weyl algebra^° of operators, W {s'^)W {s'"^) = exp[— Z2a(s='=, s'='=)]T^(s='' + s'"^). Here the 
exponent in the phase-factor q(s^, s'^) is given by 

e±e7(s±) e'±e7(s'±) 



^The zero sum condition on the matter charges stems from technicahties related to the scalar field zero 
mode [78]. 

^°The definition of the Weyl algebra follows in the standard way from the Poisson brackets between 
Y^{x),Y^{y) and an application of the Baker-Campbell-Hausdorff lemma. 
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Here a(e±,e'±) = {K,>±{f{e^)) - K,,±{h{e^))) - (ac,± (/(e'±)) - K,±{b{e'^))). Here /(e), 
6(e) are the final and initial points of the edge e respectively. Kg is defined as 

I 1 if x is in the interior of e, 
Ke(x) = < 

I 2 if a; is a boundary point of e. 

• Representation: W{s^)W{s'^) = e^^{=^a{s^ ,s'^))W{s^ + s'^). 
4.2.3 Kinematic Hilbert space 

The kinematic Hilbert space ?^kin is the product of the plus and minus sectors, ^^j^, each of 
which is a product of the appropriate embedding and matter sectors, ^^in spanned by an 
orthornormal basis of charge network states. 

A charge network state in is denoted by 18=*=) := Tg± ® W{s'^). 

The label is specified by s,± := {7(3^), {k\, . . . , {k\ , l\ )}. Here we have used 

e, e, 6,6, 

the equivalence of charge networks to set 7(3 ) := 7(5^) = 7(5'^) so that each edge of the 
charge network is labeled by an embedding charge and a matter charge. 

We refer to the quantum theory as Polymer Parametrized field theory (PPFT). 



4.2.4 Unitary representation of gauge transformations 

Finite gauge transformations generated by the density 2 constraints act, essentially, as two 
independent diffeomorphisms of the spatial manifold, one which acts only on the '+' fields and 
one which acts only on the '— ' fields. Consequently, in analogy to spatial diffeomorphisms in 
LQG, their action on charge networks is to appropriately 'drag' them around the circle^^. 
Then a unitary representation of the gauge group is given by 

?7±(0±)r,± := V(,±)?7^(0^)r,± := r,±. 

Denoting Tg± (8> W{s'^) by |s^), the above equations can be written in a compact form as 
|sj±):= ^±(<^±)|s±). 

4.2.5 Physical Hilbert space 

Physical states are obtained by group averaging the action of the finite gauge transformations 
discussed in the previous section. Henceforth we restrict attention to a physically relevant 
superselected sector of the physical Hilbert space. This sector is obtained by group averaging 
a superselected subspace, Vss of 7^ kin, T^ss = T^ts ® ^7s- 

Vfg is defined as follows. Fix a pair of graphs 7^ with A edges. Place the embedding charges 
such that k\ — k\ = for all Consider the set of all charge-network states 

{|s^) = |7^,A;^, . . . where l\ G Ze are allowed to take all possible values subject 

to the zero sum condition Yli = 0- Let Vf^ be finite span of charge network states of the 
type{|s±±^±),V(/>±}. 

^^Due to the quasi periodic nature of it is more appropriate to think of these diffeomorphisms as being 
periodic diffeormorphisms of the real hue. Consequently the action of these gauge transformations in quantum 
theory also keeps track of 'factors of 27r' when embedding charge edges 'go past x — 2n'. We refrain from going 
into such technicalities in this review. 
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The action of the group averaging map rj^ on a charge network state in yields the 
distribution 

s±e[s±] <^±eDifff+,R 

Here [s^] is the equivalence class defined by [s^] := {s'^|s'^ = for some i;^'^}, and Diff^±jR 
is a set of gauge transformations such that for each s'^ G [s^] there is precisely one gauge 
transformation in the set which maps s=t to s'=t. The space of such gauge invariant distributions 
comes equipped with the inner product 



(^±(|sf)),r/±(|sf)))phys = r?^(|sf))[|s±)], 

phys- 



which can be used to complete rj^iVf^) to the Hilbert space T-L^f^ . We shall restrict attention 



to ^phys ^phys ® "^phys- 

4.3 The construction of Tidiff 

Spatial diffeormorphism invariant distributions are constructed by the action of the group av- 
eraging map, r/(jiff on the dense space of finite linear combinations of charge network states as 
follows. Let [s"'",s~] be the orbit of s"^, s~ under all spatial diffeomorphisms so that [s"'",s~] is 
the set of all distinct charge network labels obtained by the action of spatial diffeomorphisms 
on s^, s^. Then 

?ydiff(|s+) |s-)) = r/[,+ ,,-] Yl (s'+|0(s'-|. (4.3) 

s'+,s'-e[s+,s-] 

Here ?y[s+,s-] is a constant which depends only on the orbit of s+, s^.^^ 

4.4 The Hamiltonian constraint operator in PPFT 

This section is devoted to the construction of the Hamiltonian constraint as an operator on 
the space of diffeomorphism invariant states. We follow the strategy used in LQG. Our aim is 
to first define a discrete approximant to the Hamiltonian constraint on a triangulation of the 
spatial manifold, promote the expression to an operator on the kinematic Hilbert space and 
then show that its dual action on diffeomorphism invariant distributions, admits a well defined 
continuum limit. 

Quantization of Hamiltonian constraints suffers from similar type of regularization dilem- 
mas as the Hamiltonian constraint of LQG. We would like to see if there exists any choice of 
regularization, such that 

(a) The quantum operators carry a faithful representation of the Dirac algebra (or at least 
the Witt algebra). 

(b) The kernel of the quantum constraints (or more precisely the kernel of the Hamilto- 
nian constraint plus group averaging over spatial diffeomorphisms) contains (as a vector 
space) Tiphys- 

From equation (4.1) the smeared Hamiltonian constraint with lapse N(x) is 



Ch.m[N] = j N{x) Ii+{x)X+'{x) - Ii.{x)X-\x) + + /") 



1 



yfX+'{x)X-'{x) 



^^The arbitrariness in the choice of this constant can be reduced by requiring that rydiff commute with all 
diffeomorphism invariant observables. 
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On a triangulation T, a discrete approximant to the above expression is given by 



AeT 



n+(6(A)) 



n_(6(A)) 



X+(m(A)) - X+(m(A - 1) + L5ua),o) 



lAI 



X-(m(A)) - X-(m(A - 1) - L5,(a),o) 



lAI 



+ i(y+)2(6(A)) + i(y-)2(6(A)) 



(4.4) 



where b{A) is the beginning vertex of simplex A, | A| is its length and m(A) its midpoint. The 
symbol A — 1 denotes the simplex to the left of A, and it is understood that if Ai is the left-most 
simplex with b{Ai) = then m(Ai — 1) = m{Ai\f), with A^v being the right-most simplex such 
that /(Ajv) = 27r where /(A) is the ending vertex of A. 

Since only the holonomies of II-i-, are well defined operators on Hkin, the local fields II-i-, 
^y±^2 j^gg^ |-,g approximated on T by appropriate combinations of holonomies. The equa- 
tion (4.4) is a discrete approximant to the continuum expression and we may modify it by terms 
which vanish in the continuum limit, |A| — t- 0. It is straightforward to see that the following 
expression is one such modification, 



-ihN{b{A)) 



AI2 



1 + 



ih 



-n_(6(A))(6(A)) 



l + ^n+(6(A))X+'(6(A)) 



-ih 



l + ^iY+f{b{A)) 



-Aih' 



1 + 



. (4,5) 



In a nutshell our aim is to choose a quantization scheme for Cham.T [N] such that the quantum 
constraints are a faithful representation of the classical constraint algebra, and the kernel of 
Hamiltonian constraint (after averaging over spatial diffeomorphisms) leads to ^phys- One 
way to anticipate what kind of regularization choices one should make is to consider action 
of classical Hamiltonian vector field of the hamiltonian constraint on cylindrical functions. If 
we start with the density one constraint, this action is quite involved due to the presence 
of factor. However as we will see below, in quantum theory the spin-network states are 
eigenstates of inverse metric operator, and the non-triviality of quantum action basically comes 
from n+X+ — n_X~ + (y^)^ + (y~)^. whence let us focus on that part, which is really just 
the density two constraint. 

It is a simple computation to show that the action of Hamiltonian vector field of (density 
two) constraints on cylindrical functions can be written as 



= lim \ {Tg+ \ 



)*H+]-r,+ [n+]), 
)*n_]-r,+ [n^]). 



Whence we seek finite triangulation holonomy approximants to the various local fields of 
interest in such a way that Cham,r[-^] is proportional to a combination of finite gauge transfor- 
mations minus the identity, with the finite gauge transformations being parametrized by |A| so 
that at I A| =0 the gauge transformations are just identity. 
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We display exactly such approximants to II-i-, (Y^)"^ in next section. The holonomies turn 
out to be state dependent^'^, the approximants yield n-i-(6(A)) and (y^)^(6(A)) to leading order 
in |A| and the resulting quantum constraint is a linear combination of operators localised on 
vertices, each of which is proportional to a the difference of a finite gauge transformation and 
identity. 



4.5 Suitable "curvature approximants" 



We shall focus on the left moving (+) variables. Given a state Tg+, let Ve{'j{s~^)) be the set of 
all vertices v such that A:^ — /cj, 7^ (where e„(e^) refer to edges which terminate (originate) 
at V. Similarly let Vivf(7(s^)) be the set of all vertices v such that — / 0. Then as 
shown in [77], the quantization of n+(w) which defines the Hamiltonian constraint whose kernel 
contains states in "Hphys is 



n+(^)T,- 



|A|(A;+ 
i 

\A\ 



{hi - m- 



\iveVE{l{s+)), 



(4.6) 



Next, we show that the above choice (4.6) directly leads to an operator action of the '+' embed- 
ding part of the Cham,T which is a finite diffeomorphism on the + part of the embedding state. 
As we shall see, this will finally lead to a satisfactory definition of the Hamiltonian constraint. 
With the above definition, the approximant to the II+X^' term is 

n+X+'(6(A))|Tr,+ := ^flJ(6(A))(X+(m(A)) -X+(m(A- l)))r,+ 



A 



lAP 



i)r,- 



if 6(A) G Ve{i+ 
if 6(A) ^ Ve{i' 



Whence for x E Ve(7^ 

n+x+'(6(A))|Tr,+ 



|A|2 
—ih 



AP^^4, 



-ih 



lAP 



(t/+'^(<AA) - l)Ts- 



Here, (j)/^ is a diffeomorphism of the circle (more precisely, (j)/^ is a periodic diffeomorphism of 
the real line) which is identity in the neighborhood of all the vertices of T except 6(A), /(A).^^ 
Further, maps 6(A) to /(A) and its action on the charge network label s"^ is denoted 
by s~^^. U~^'^{4'a) is the restriction of the unitary action of the finite gauge transformation 
[/+(0+ = to the left-moving embedding Hilbert space. 

Similarly the quantization of (y+)2) that yields a satisfactory definition of the Hamiltonian 
constraint is 



{Y+nb{A))TW,- 



exp 



4k 



/ I+l 



"A 



^^As we will see shortly, this dependence involves the eigenvalues of the embedding operators X^. Since 
are the analogs of the LQG densitized triad operators, this feature is reminiscent of the p, scheme employed in 
the improved dynamics in LQC [9]. 

^"^In addition, (j>/s also differs from identity at x 
a consequence of the circular topology of space. 



2tv (or a; = 0) if 6(A) = (or /(A) = 2n)- this is just 
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which can also be rewritten as 

-m 



(?+r^(&(A))T=^(f>+'^^(0A)-i) 



where has been defined above and U~^'^ {(t>/\) is the restriction of the finite gauge transforma- 
tion operator lJ^{(j)'^ = to the matter Hilbert space. The analysis of the right-moving (— ) 
mode proceeds in a similar way. 

Substituting the above "curvature approximants" into (4.5) leads to the following definition 

of Cha.ra,T[N] On Tikin, 

Cham,r[A^]|s+,s-) = Yl N{b{A)) 

AeT 

b{A)&VE{s+)UVE{s-) 

X "?7+'^(0a) ® U-^^{<Pa' i) ® U+^'H^a) ^ U-^^'i^A i) - ll |s+,s-). 

^ ^ \A\y/X+'X-' 

Using the fact that the unitary operators in the above equation are just restricted actions of 
unitary operators associated with finite gauge transformations, we obtain 

Cham,T[iV]|s+,s-) = N{b{A)) A(s+, , 6(A)) 



AeT 

b(A)eVB(s+)UVi;(s-) 



\s , s 



This action precisely mimics the action of the Hamiltonian vector field of Cham[-^]- We could 
now consider the continuum limit of Cham.T [N] on "Hkin in the URST topology. As shown in [77] , 
just like in LQG, this operator converges to a continuum operator ChamlA^] on T^kin- 

Specifically, in the notation used above, we may choose the limit of the one-parameter family 
,T{'y,5)W] to be the operator Cham[A^] where 

C'ham[A^]|s+,S~) := C'ham,T(7A,)[^]|s^'S~)- 

On one hand this is a very promising candidate for the Hamiltonian constraint of PFT. A moment 
of thought reveals that given any (^'| in Tiphys 

(^'|C'ham[iV])|s+,S-) =0, 

for all and Is"*", s^) . On the other hand, this Hamiltonian exhibits the same triviality properties 
that the LQG Hamiltonian does. Instead of going into technical details, we summarise the 
findings of [77] in qualitative terms. 



4.6 Commutator of Hamiltonian constraints in PPFT 

(1) Precisely due to the reasons we elaborated upon in the earlier sections, the commutator of 
density one constraints C'ham[A^] vanishes on T^kin- 

(2) One can define an analog of the LM habitat, on which both the r.h.s. and the l.h.s. 
of constraint algebra can be quantized. Due to the density one character of the constraints, 
the product of the two Hamiltonians is ultra-local in the lapses and hence vanishes under anti- 
symmetrization. It turns out that if we use the same curvature approximants that we used to 
define Cham,T[A^] to define r.h.s., then two things happen. On one hand the operator at finite 
triangulation Cdi{{,T[Q~^[^^ ^]] mimics the action of corresponding classical Hamiltonian vector 
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field, and on the other hand Cdiff,T[9~^[-^) -^]] is proportional to |A|, whence on LM habitat, 
its continuum limit is zero. 

(3) As we had noted earlier, in their analysis of the LM habitat in LQG, the authors of [64] 
noted certain necessary conditions for the commutator to be non-trivial. One of them was the 
density 7^ 1 character of the constraint, and the other was, a definition of curvature approximant 
such that the newly created vertices where themselves not in the "kernel" of the Hamiltonian. 
The choice of curvature approximants in PFT ensure the second condition is satisfied. So one 
would hope that appropriately changing the density weight (which is dictated by how many 
powers of | A| we need in the denominator, so that the continuum limit contains correct number 
of derivatives) would give us a non-trivial constraint algebra. 

In two dimensions, in order to get the right powers of |A| in the denominator, one needs to 
work with density two constraints. But when we work with the density two constraints, several 
remarkable things happen. 

(3i) The continuum limit of Cham,T[-^] (we use tilde to denote the density two character) is 
not well defined on ^kin or the LM habitat. However, the commutator of two such Hamiltonians 
is well defined on the LM habitat and has an anomaly! 

{3ii) There exists a different habitat, on which the density one constraints admit a continuum 
limit and the constraint algebra closes off-shell without any anomaly! and all the states in "Hphys 
are contained in this new habitat. 

Whence we see that at least in the case of our toy model, everything nicely fits into a coherent 
picture. The requirement that constraints close off-shell is a very stringent requirement. It not 
only pins down the curvature approximants and the density weights, but also dictates that the 
right home for the constraints might be something quite different then 'Hkin or even Tidm- In 
the case of PPFT, it turns out that the kernel of such constraints rather remarkably contains 
the solution obtained via completely different methods. 

5 Diffeomorphism constraint in LQG 

If the lessons learnt above in the case of PPFT could be applied to LQG, one could in principle 
obtain a quantization of the Dirac algebra. The key lessons one learnt in PPFT were 

1. Higher density objects could be well defined in the theory if one worked on suitable habi- 
tats. 

2. In the case that the constraints generate a true Lie algebra, the constraint operators at 
finite triangulation (with triangulation parametrized by 5) should be related to the action 
of finite gauge transformation if the finite gauge transformations were unitarily represented 
in the theory. 

It is unclear at this stage whether these lessons could be applied to the Hamiltonian constraint 
but as a warm-up one could apply these to the diffeomorphism constraint. 

Recall that in the Dirac algebra, there is a Lie subalgebra generated by spatial diffeomor- 
phisms, 

{HdiH[N],HdiH[M]} = -Hdm[N,M]. (5.1) 

We can now ask if there exists a quantization of //diff[^] on some subspace of Cyl* such that 

(a) The operator Hdm[N] is a faithful representation of (5.1) on that subspace. 

(6) The kernel of -ffdiff[-^] matches precisely (or at least contains) the vector space of diffeo- 
morphism invariant states obtained via Rigging map techniques. 
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Figure 2. 

It is important to note that at the level of finite triangulation, HTdmi^] is well defined 
on 'Hkin and also admits infinite number of quantization choices (choice of triangulation, choice 
of loops to define operator for curvature, choice of surfaces to define operator for triads, choice 
of representation for the curvature operator, etc). Whence it is interesting to see if requiring the 
diffeomorphism constraint operator to satisfy (a) and (6) restricts these choices in anyway. Also 
as the r.h.s. of Poisson bracket of two Hamiltonian constraints has a diffeomorphism constraint, 
it might be essential to have an operator correspondent of H^mflN]. 

As shown in [79], there does exist a quantization of -ffdiff[-^] which precisely generates 
infinitesimal diffeomorphisms on LM habitat, and hence satisfies (5.1)^^. 

The ideas underlying the quantization are sketched below. As always we set G = h = c = 1. 
We shall also set the Barbero-Immirizi parameter to unity. The diffeomorphism constraint 
Hdm{N) is 

H^m{N)= I C^A^,Et = V{N)-G{N-A^^), 
Jt, 

where V{N) = N^F'^^E^ and g{N^Al) = N^AlVaEf. 

Let T{6) be a one-parameter family of triangulations of S with the continuum limit being 
(5 — )• and let -ffdiff ,T(<5) ) ^T(5)i Gt{S) be finite triangulation approximants to the quantities -ffdiff; 
V, Q of the above equations. Thus Dx^s)-, ^T{S)-, Qt{S) a-re expressions which yield /^diffj V ■, Q in 
the continuum limit. 

For simplicity consider a (non- gauge invariant) spin network state consisting of a single 
edge e with spin label j so that the state is just the m, n component of an edge holonomy hem" 
of the (generalized) connection along the edge e in the representation j, the indices m, n taking 
values in the set 1, . . . , 2j + 1. In what follows we shall suppress some of these labels and denote 
the state simply by hg. 

As we want the operator at finite triangulation to equal finite diffeomorphism minus the 
identity, our first goal is 

(l + i5HdiH,T{5))he = ^0(5,jv)oe' 

where (p{5, N) o e is the image of e by the diffeomorphism <j){5, N) which translates e by an 
amount 5 along the integral curves of the shift vector field N'^ (see Fig. 2). 
As shown in [79] the desired result is obtained in the following steps. 

(i) First we set 

(1 + i(5Fdiff,T(5)) := (1 - »^^T(<5))(1 + i^ynS))- 

^^Recall that on Hmu only finite diffeomorphisms are unitarily represented and the generator of such 
diffeomorphisms does not exist. However such generators could and indeed they do exist on certain distribu- 
tional spaces. 
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Figure 4. 



(ii) Next, one shows that 

{l + i5VT(5))he = h-^j^^^y (5.2) 

Here e(A^, 6) has the same end points as e (as it must by virtue of the gauge invariance 
of V) and is obtained by joining the end points of (p{5, N) o e to those of e by a pair of 
segments which are ahgned with integral curves of A^'* as shown in Fig. 3. 

(iii) Finally, one shows that the Gauss Law term, (1— removes these two extra segments 
(see Fig. 4). 

The major part of the analysis concerns the derivation of the identity (5.2) in step (ii) above. 
V'p(^S) is written as a sum over contributions Va where A denotes a 3-cell of the triangulation 
dual to T((5), and is a finite triangulation approximant to the integral N"'F^f^E^. 

Choice of operator ordering. Order the triad operator to the right in Va so that only 
those 3-cells contribute which intersect e. 

Choice of triangulation. The triangulation T{5) is adapted to the edge e so that its 
restriction to e defines a triangulation of e. Thus, there is a triangulation of e by 1-cells and 
vertices of T{6) so that each of these vertices vj, I = 1, . . . , is located at the center of some 
3-cell A = A/. 

Choice of curvature approximant. We define a finite triangulation approximant to F^^ 
in a such a way that the following identity holds: 

(l + idVAj)he = /lg(A_r) 

V/ G {1, . . . , — 1}. Here e(A/) is obtained by moving the segment of e between vi and 
along the integral curves of N"" by an amount 6 and joining this segment to the rest of e at the 
points vi, by a pair of segments which run along the integral curves of A^" as shown in 
Fig. 5. 

This can be accomplished by using the full quantization ambiguity that is available to us 
when defining curvature operator in LQG. More precisely not only does one have to tune the 
spin (in which the trace of holonomy is evaluated) to the spin of the underlying state, but one 
also has to add the "higher-order operators" obtained by involving the fluxes. For more details 
we refer the reader to [79]. 

Sum to product reformulation: Finally the contributions from all the V/\j yield the edge 
e{N,6) of Fig. 3. Recall that Vx(s) is obtained by summing over all the cell contributions Va. 
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Figure 6. The double lines indicate retraced paths. 



However, summing over the action of all the Va^ on he only yields a sum over states of the 
type /le(Aj)- III order to obtain the desired result, h-^j^ the sum over A is first converted to 
a product over A, i.e. to leading order in S, we show that 

l + iS^VA^Ylil + idVA). 

A A 

Hence, replacing the sum over the corresponding operators by the product provides an equally 
legitimate definition of Vj'(s)- The replacement then leads, to the following identity 

N-l 

J](l + i6VA)he = n (1 + ^'^^A,)/ie. 
A 1=1 

The definition of Vaj is such that each factor in the product acts independently, the Ith factor 
acting only on the part of e between vj and We are then able to show that the result (5.2) 

follows essentially through the mechanism which is illustrated schematically in Fig. 6. 

Thus we see that by making some very specific choices of triangulation, curvature and ope- 
rator ordering, one obtains 

^diff,T(5)[^]|7, = -^{U{<pf) - l)h,lS). 

It can be easily shown that -f^diff ,T((5) [^] admits a continuum limit on LM habitat where it 
precisely generates Lie action by N. 

6 New progress and insights from the BF topological model 

We now present a different family of toy-models, namely the topological BF field theories in 
arbitrary dimensions. It is a field theory which has the same unreduced phase space as Yang- 
Mills or general relativity with Ashtekar-Barbero variables. But as a topological theory it is 
an exactly solvable theory and it has been considered as the main testing ground for the loop 
quantum gravity program since its early days, and even more after the birth of the spin foam 
program. Just like in the previous sections, the idea is to quantize d la loop the unreduced 
phase space, and then impose the topological symmetries at the quantum level. One expects 
to reproduce in this simplified situation a quantum algebra of constraints, and the resulting 
quantization should be equivalent to the one which has been performed using standard path 
integral methods. 
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This section first reviews the BF theory itself, in its standard quantization, then the most 
recent works in the loop framework, emphasizing the geometric aspects which are expected to 
show up similarly in four-dimensional quantum gravity. A nice review from the loop and spin 
foam perspective is [10] . The present section should be considered as a new review on the state 
of the art for the loop approach to this topological model, ten years after. 

We focus on the symmetry group SU(2), but arbitrary spacetime dimension n. The BF 
field theory was introduced by Blau and Thompson in [25] and by Horowitz in [72]. It also 
appears in a number of interesting field theories. When n = 2 [114] and n = 4 [40], it appears 
as the zero coupling limit of Yang-Mills theory, an approach known as BF Yang-Mills (and 
as such it was observed early in lattice gauge theory). In three dimensions, it describes pure 
Riemannian gravity (with degenerate metrics) , and remarkably can be recast in a Chern-Simons 
formulation [113]. In four dimensions, it is also possible to establish a link with gravity (which 
is crucial to spin foam models for quantum gravity) through the BF Plebanski formalism [43], 
or the BF MacDowell-Mansouri formulation [62]. 

It is also an interesting model exhibiting a relationship between field theory and topological 
invariants. The partition function produces the Ray-Singer torsion [24], and expectation values 
of observables are studied in [39]. 

6.1 Transition amplitudes 

6.1.1 Field equations and flat connections 

Let P be a principal SU(2) bundle, over a n-dimensional manifold M. We consider a connec- 
tion A on P, locally seen as a su(2)-valued 1-form over M with specific transformation rules. 
We also introduce a (n — 2)-form taking values in the adjoint bundle, i.e. B is su(2)-valued 
and transform under the adjoint action of SU(2). When d = 3, B physically corresponds to the 
(Riemannian) triad, usually denoted e, and the spacetime metric is obtained from g^i, = Jjje^ei. 

To write the BF action, we consider the non-degenerated bilinear invariant form 'tr' on the 
algebra, and 



with Fa = dA + g [A, A] being the curvature of A. The group G of usual gauge transformations 
consists in G-valued fields, which generate the following transformations, 

A^ AdgA + gdg~^, B^AdgB. 

In matrix notation the adjoint action reads Adg X = gXg^^. Since 5bf is the integral of 
a differential form, it is naturally invariant under diffeomorphisms of M. 
The equations of motion are particularly simple, 

dAB = Q, Fa = 0, 

where dA = d + [A, ■] is the covariant derivative. It turns out that all solutions are locally 
identical up to gauge transformations. Notice that the second equations comes from the fact 
that B is a. Lagrange multiplier for the curvature, so that the connection has to be flat. We 
know then that all flat connections are locally gauge equivalent. We also get that B is locally 
exact, B = dA^l^- Now it is time to state that the action has an extra invariance under 

A^ A, B B + dAV, 

for any (n — 3)-form r] over M taking values in su(2). That additional gauge symmetry is often 
called translation symmetry. Denoting T the group of such translations, the full group of gauge 
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symmetries is the semi-direct product Q <k T. It is then possible to transform any solution B 
into the trivial solution. The physical consequence is that the theory is deprived of local degrees 
of freedom. 

Those equations admit a useful cohomological interpretation. Squaring the covariant deriva- 
tive gives d\oj = [-Fyij'^] which vanishes on-shell. So dA defines cohomology classes. The set 
of non-gauge-equivalent solutions to the classical equations of motion consists in non-gauge- 
equivalent flat connections with elements of the (n — 2)th cohomology group, the latter being 
isomorphic to the second cohomology group. 

It is interesting to note that general relativity and the BF model share some essential prop- 
erties which enable to make geometry dynamical. They are both formulated using differential 
forms, without any spacetime metric, and are diffeomorphism invariant. Hence, it should be 
clearly understood that what makes the BF model special is the translation symmetry, so that 
diffeomorphisms are actually encoded into the full set Q kT. Suppose ^ is a vector field on M, 
then using the inner product zg, we can form the function (p^ = i^A and the (n — 3)-form 
r/^ = i^B. The translation symmetry and diffeomorphisms are related in the following way 



where is the Lie derivative along ^. When written in this way, infinitesimal diffeomorphisms 
appear as gauge transformations up to terms proportional to the variations of the action 
6S/5B = F, so that the two symmetries are de facto equivalent. Obviously in general rela- 
tivity, it becomes impossible to rewrite diffeomorphisms as gauge symmetries plus terms which 
vanish on-shell. 

In the context of the BRST (or Batalin-Vilkovisky) formalism, the symmetry of the BF 
model has been interpreted as a vector supersymmetry [65, 86], which turns out to be crucial 
to solve the local cohomology of the BRST operator, and in the end to prove the perturbative 
renormalizability [85] ! An important aspect in that algebraic approach is that the BRST algebra 
does not close off-shell. The BRST operator is nilpotent only on-shell when n > 4, due to the 
fact that the translation symmetry is on-shell reducible. Indeed, the transformation of B is the 
same for r/ and r} + dA'fp- As a consequence, care is due when trying to quantize it in dimension 
four and higher. 

6.1.2 The Hamiltonian analysis and flat connections again 

The Hamiltonian analysis is straightforward. One finds that one canonical variable is the pull- 
back of the connection to the canonical surface S, denoted A\. The other canonical variable is 
the analog to the electric field, i??, transforming in the adjoint representation. Those variables 
satisfy the canonical bracket 



The phase space is actually identified as the cotangent space T*A to the space of connections A, 
exactly like in Yang-Mills theory. 

However, the Hamiltonian is different since it is a combination of constraints imposed with 
Lagrange multipliers. The constraints are simply the restrictions of the equations of motion 
to S, which do not involve any time derivatives. We integrate the local constraints as follows 



5(^A = Li^A = dA(l)( + i^F, 



6^B = C^B = dAr]^ + [05, B] + ^ dAB, 



{E^{x),Ai{y)}=6^6l6(^Hx-y). 




■)) 




where D is the covariant derivative with respect to j4^. Through the Poisson brackets, those 
constraints generate the whole set of gauge transformations on the canonical variables. Hence 
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they satisfy the expected algebra, 

{Cg(A),Cg(A')} =Cg([A,A']), {CG{A),CFiN)}=CF{[A,N]), 
{Cf{N),Cf{N')}=0, 

which is recognized as the Poincare algebra ISU(2). The algebra is first class and the local 
number of degrees of freedom is easily evaluated to zero, recovering the result of the cova- 
riant analysis. Cq is the Gaufi law which generates local SU(2) rotations, while Cp generates 
translations as expected since curvature is zero. 

To completely solve the theory classically, one can first reduce the phase space by the Gaufi 
law Cg, just requiring invariance under the standard SU(2) gauge transformations. That leads 
to the reduced Yang-Mills phase space T*{A/G). Reducing with respect to the constraint Cp 
implies that we consider only the set of flat connections. It is also necessary to identify 
any two momenta which differ by a transformation generated by Cp, i.e. Ef ~ i?? + e"'^'^ D^-qc. 
The relevant set of momenta is thus the second cohomology group of D (which squares to zero 
on ^o)- In conclusion, the reduced phase space is the set of flat connections modulo the ^-action 
together with the second cohomology group over each flat connection. This is the set T* (Aq/G), 
as described in [72, 115]. 

There exists a description of the set of flat connections which shows that the reduced phase 
space is often finite dimensional and which proves useful in practice. Holonomies of a flat 
connection only depend on the homotopy types of its paths. Hence a flat connection defines 
a morphism from the fundamental group vri($]) into the structure group SU(2). Under a gauge 
transformation, such a morphism gets conjugated by the value of the transformation at the 
reference point of 7ri(S). We thus get 

Ao/g C Hom(^i(S),G)/G, 

where G acts by conjugation. The set Hom(7ri(I]), is known as the moduli space of flat 

SU(2)-bundles. When 7ri(S) has a finite presentation, then it is a real algebraic variety. It also 
usually has singularities which prevent it from being a smooth manifold. Later, we will ignore 
those singularities and consider the case where Ao/G really coincides with the moduli space of 
flat bundles. 

This space has been mainly studied when S is a two-dimensional closed surface, where a sym- 
plectic structure is available. That has been quite efficient to solve Yang-Mills theory in 2d [114], 
Chern-Simons theory in the canonical formalism [113]. More recently [22, 35, 36, 63]. 

6.1.3 Quantization 

The path integral. Formally, without gauge fixing, it reads 

"Z(M) = / DADBexpU [ tr(5 AF^)) = [ DA TJ '^(^^(a;))", 

from which we learn two things. 

• The functional integral is exactly peaked on flat connections over M. 

• Gauge fixing is crucial to get a topological invariant (if any). Indeed, if we naively continue 
the above calculation, we find (after some regularization) the determinant of 6Fji/5A which 
is the the covariant derivative acting on 1-forms and which is not a topological invariant. 
The intuitive reason is that dA here only probes the topology going from 1-forms to 2- 
forms. The final topological invariant should instead probe all dimensional forms (or cells) 
over M. 
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The path integral can be evaluated using the usual Faddeev-Popov gauge-fixing in three 
dimensions, and the Schwarz' resolvent in higher dimensions (to deal with the reducibility), 
see [24, 25]. 

Let us sketch the three-dimensional situation. We assume for simplicity that the moduli 
space of flat connections is finite and that for any flat connection A(^a) the cohomology of the 
operator ^Aj^j is trivial. We pick up a background metric g and select a representative of 
each equivalence class of flat connections up to gauge which satisfies c?A(q) * ^{a) = 0; where * 
is the Hodge operator. Setting A = A^^^-^ + a in the neighborhood of we get 

Z(M) = ^/x(^(„)), 

a 

with 

^{A^^-j) = j DaDBDuDODxDxDuDv 

X exp^i j tr (5 A dA^^^a + u A *a + v A dA^^^ * B + oj A Aa^^^uj + x A Aa(„)X)^ , 

where A^^^j = dAj^^^ * dAf^^y u, v are Lagrange multipliers imposing gauge conditions on a, B, 
and uj, X O, x a-re the corresponding ghosts and anti-ghosts. We refer to [24, 25] for details. 
The main point is that ^(^(q)) can be evaluated to be the Ray-Singer analytic torsion. It is 
a topological invariant, which does not depend on the choice of the background metric. 

The wave-functions. From the canonical analysis on M = S x [0, 1], it is clear that they are 
square-integrable functions over the set of flat connections modulo gauge transformations [72]. 

The measure. The inner product is obtained by considering the path integral on M = 
S X [0, 1]. The measure is given by the Ray-Singer torsion on M which in this case descents 
on S, so that the measure is well defined on the canonical surface. 

A physical IR divergence. The partition function is well-defined for topologies where the 
cohomology H*(dA) is trivial. At least in three dimensions, there is a nice physical interpre- 
tation, noticed by Witten [115]. There the partition function diverges after gauge fixing if the 
moduli space of flat connections is not a finite set. Notice that the equation for the tangent 
vectors to that space is SF = dA^A = 0. This is exactly the equation of motion for the triad 
dAG = 0. Hence, the partition function is finite on those topologies where there are no non-trivial 
(equivalence classes of) solutions for the triad. The intuitive picture of the quantum mechanical 
consequence is that there are only some quantum fluctuations around the trivial solution. When 
restoring the units, it means that lengths are fluctuating around the Planck scale, explaining 
the finiteness of Z{M). However, if there are some tangent directions to the moduli space of 
flat connections, there are some non-trivial (equivalence classes of) triads satisfying d^e = 0. 
The integrals along those directions are non-compact, hence the divergence. Physically, that 
corresponds to an IR divergence, since the main contributions to the path integral come from 
large triads, hence spacetimes much larger than the Planck scale. The result is that a macro- 
scopic physics can exist only on topologies with non-trivial cohomology and a divergent partition 
function. 

6.1.4 Loop quantization 

Since the loop quantization relies on cylindrical functions supported on graphs and since the 
measure on the set of flat connections is the torsion, it would be convenient to evaluate the torsion 
using cell decomposition. That notion exists and is known as the Franz-Reidemeister torsion. 
This actually shows that the BF model admits an exact quantum formulation in a discrete 
setting. Some details can be found in [37]. 
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Discrete flat connections and covariant derivative. The main idea is to find an 
equivalent on cells to the covariant derivative. Let us describe its parts which act on 0-cells and 
1-cells. Further details can be found in [22, 35, 36], in [63] in the context of knot invariants, and 
also in the Chern-Simons literature (where the measure is rather the square-root of the torsion) 
such as [73, 101]. 

Let A be a cell decomposition of S, and Aj the set of i-cells. Since we focus on holonomies, 
we can forget that they come from a connection, and consider instead the notion of discrete 
connection which is the assignment of SU(2) elements to edges of A, denoted A = (s'e)eeAi- 

Gauge transformations induce a SU(2)l'^ol-action. A discrete gauge transformation h is there- 
fore a set of group elements (/i^)^eAo acting at the vertices of A according to 

^•^= (V)9e^;(e))eeAr 

To evaluate the effect of a small gauge transformation, we consider the differential at /i = 1 of 
the map which sends h to h ■ A, 



d^A\i - su(2)l^oUrASU(2)l^il, {Rg^, {id- Adg^,)v) 



eeAi 



where Ad stands for the adjoint representation of the group on its algebra, and for the 
right translation. In matrix notation: AdgV = gvg~^ and Rg^,v = vg. The kernel of d'^A\i is 
the algebra of the isotropy group C(^) of the connection A, while its image corresponds to the 
directions along which A is changed by the group action. If Oa denotes the orbit through A, 
then we know that 

Oa^SU(2)/C(A), and TAOA = i^d^A\i- 

The curvature of a connection is encoded in the holonomies along the boundary of 2-cells. It 
is represented by the map 



H : SU(2)l^il ^ SU(2)I^2|^ ^ ^ (^/(^) = JJ 5^^') 



where eje is the relative orientation of the face / on the edge e, and Hf{A) is the 'holonomy' of 
the connection A around the face /. This provides a notion of flatness on A: the connection is 
flat if 

H{A) = 1. 

For a discrete flat connection cp, the tangent space to the moduli space of flat SU(2)-bundles 
hom(7ri(S),SU(2)) is 

r^(hom(7ri(S),SU(2))/SU(2)) = ker dH^/ imd^^\i. 

Moreover, gauge invariance implies that dH^od^y^^i = 0. That equation means that cohomology 
classes can be defined. The analogy with the continuum is that d'y^\i is the covariant derivative 
on su(2)-valued 0-form, and dH^ the covariant derivative on su(2)-valued 1-form. Hence, the 
tangent space hom(7ri(S), SU(2))/ SU(2) is the first cohomology group. That provides with 
an algebraic description of the local properties of the moduli space of flat connections. 

Reidemeister torsion. The above description gives a (co)chain complex which includes 0-, 
1- and 2-cells. If A is a cell decomposition of dimension d greater than two, the complex can be 
extended for any discrete flat connection (/>, 

■J^^^CA) ^ <5i(A) (50(A) 

^ Q^(A) ^-^^ ...< C^(A) C^(A) C°(A) ^ 0, (6.1) 
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where each cochain group basically contains assignments of a su(2) element to each cells, 
C^(A) ~ 5u(2)l^'l. We have also relabel the previous 0, 1-cell covariant derivatives, 5^ = d'^if,\i 
and 5]p = dH^. 

Assuming the cohomology of the above complex is trivial, one can define the Reidemister 
torsion as follows 

tor,(S) = n(det(5^t5^))(-i)V2, 

Here f denotes the adjoint with respect to the natural 5u(2) inner product. Obviously the 
operators (5^^ 5^ have some kernels so one must take care of removing them by restricting the 
operators to the orthocomplements of their kernels. Then the torsion only depends on the 
equivalence class [(/>] of the flat connection modulo gauge and is a topological invariant in the 
sense that it is independent^^ of the choice of cell decomposition A. 

Transition amplitudes. That combinatorial torsion provides a natural analog to the Ray- 
Singer torsion which is well adapted to cylindrical functions. The latter depend on a finite 
number of holonomies, just like those used to define the torsion. Moreover, provided we consider 
gauge invariant wave-functions, their value on a flat discrete connection will actually only 
depend on the equivalence class [(/)]. 

Hence, a natural candidate for BF transition amplitudes is built by gathering the torsion to 
get the measure and cylindrical functions evaluated on equivalence classes of flat connections as 
quantum states. Consider two graphs ri^2 supporting two cylindrical functions 2J1 1- Choose 
a cell decomposition of A which contains ri^2- Then the natural extension of the formal path 
integral is 

(Vt.jJVt,,/,) = 5^/i(M)tor[<^](S)/2([</.]). 
[01 

Notice that the dependence of /i,2 ([</*]) on their graphs is only through their homotopy classes. It 
implies for instance that different knotting classes of a graph embedding cannot be distinguished 
by that invariant. 

Such a formula defines the physical inner product on the reduced Hilbert space, or equiva- 
lently the matrix elements of the projector onto those states which satisfy the flatness condition 
at the quantum level. 

Obviously that proposal should be made more precise by building the full path integral on 
spacetime directly from a cell decomposition, identifying the gauge symmetries in that setting 
and gauge-fixing them appropriately. This has been done at the covariant level when spacetime 
is three-dimensional (with boundaries): the full set of gauge symmetries is described in [59] while 
the relationship to torsion can be found in [22]. At the canonical level, i.e. on a two-dimensional 
surface E, the transition amplitudes have been evaluated in [88] by defining the projector onto 
physical states which satisfy flatness. In that situation, there is no translation gauge symmetry 
to take care of, and the torsion is trivial, as well-known in Id Yang-Mills theory [114]. 

Coupling to extended matter. As a topological field theory, the BF model is known 
to feature exotic statistics [12, 23], which depending on the dimension involve particles, strings 
and/or branes. The dynamics and quantization of BF theory coupled to such extended matter 
are presented in [11, 56]. 



^^The full statement is that torsion is invariant under simple homotopy transformations of A. The notion of 
simple homotopy is somehow intermediate between homeomorphisms and homotopy. 
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6.2 The Wheeler— DeWitt equation for quantum flat space 

We showed in the previous section that the BF field theory admits an exact discretization, 
which is well-adapted to cylindrical functions, and hence to the loop quantization. If A is a cell 
decomposition of the canonical surface S, the projector on physical states, satisfying the flatness 
constraint Cj, reads formally 

Pa= U 5{Hf). 
/6A2 

Let us say a few words on that projector. When written in this way, it only sees the 2-skeleton 
of A, so that one can try to study it on a arbitrary 2-complex, which may or may not be the 2- 
skeleton of a higher dimensional cell complex. Technically, divergences arise because they are too 
many delta functions in the above expression. If only a 2-complex is given, the complex (6.1) 
is truncated to cells of dimension lower than two. Since gauge symmetries come from the 
action of dimensions three and higher onto the faces, divergences cannot be understood for the 
presence of gauge symmetries. However, the structure of the divergences are still encoded into the 
cohomology twisted by flat connections (truncated at dimension 2), and the rate of divergence 
can be extracted as exactly the dimension of the cohomology group H'^{A) (see [35, 36] and [34] 
for a summary and a comparison with other approaches to those divergences). This is an 
important result since the divergences are now well understood even when there cannot be 
interpreted as due to gauge symmetries. That may be helpful for the flat sector of full quantum 
gravity, since the gauge symmetries of BF are generically broken in the full theory. 

We now restrict to situations where this projector is well-defined. That means either the 
gauge-fixing is trivial, like when A is a decomposition of a closed orientable two-dimensional 
surface of genus g greater than two, or it is almost trivial in the sense that it is sufficient to 
remove a given number of faces which carry redundant delta functions. That is the case of the 
spherical topology, and such a gauge-fixing is described in [59]. 

The projector Pa on Tig has been constructed in [88]. Its matrix elements provide a definition 
of the physical inner product. The projector can be unfolded in spacetime, using for instance the 
spin network basis. That mimics a time evolution for spin networks which actually corresponds 
to a spin foam model, namely the Ponzano-Regge model. 

6.2.1 The new Hamiltonian for 2-(-l gravity 

Though the result of [88] was a major progress, it was still unclear whether a Wheeler-DeWitt 
equation for 2+1 gravity in the spin network basis can be written. Since the BF model describes 
in (2+l)-dimensions pure gravity (with degenerate metrics), the canonical formalism can be 
performed differently than in Section 6.1.2. It can be done like in general relativity, by intro- 
ducing the shift and lapse. They are respectively Lagrange multipliers for the vector constraint 
which generates diffeomorphisms of the canonical surface and for the scalar constraint H which 
generates time reparametrizations. 

So we can trade the constraint Cp for those two. Physical states have to satisfy the vector 
constraint Va\^) = and the Wheeler-DeWitt equation 

^IV') = 0, 

for some regularization of Va and H. There are a few things we know about the Hamiltonian, 

• from the previous section, we know we can work on a single cell decomposition of spacetime, 
because H should not change the graphs of cylindrical functions, 
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Figure 7. One gets the new Hamiltonian by contracting tlie Wilson loop around the face with the 
fluxes meeting at the starting node of g/. 

• but K should change spins in the spin network basis, so that it probably should generate 
difference equations. 

Note that on a two-dimensional surface S the curvature tensor is just a vector Fyi. 
Since i?^ is the co-triad, the vector n = Si x i?2 is the normal to the canonical surface. Then 

H = n ■ Fi2/\n\. That is the origin of most difficulties to regularize H on cylindrical functions, 
as we explain. 

Assume that A is a complex dual to a triangulation. The natural variables are the Wilson 
lines along edges of A and their canonical momenta El which are fluxes through those edges. 
The geometric picture is that the fluxes represent the edges of each triangle of the triangulation 
embedded into M^. The curvature at a node of the triangulation is regularized by the Wilson 
loop around the face dual to the node. However, there no clear definition of the normal n to 
the node. 

The way-out which has been proposed in [29] is the following. If the node of the triangula- 
tion is for instance 3-valent, it is shared by three triangles, and one can define three regularized 
quantities looking like H using the normal to each triangle. Since the three normals are generi- 
cally linearly independent, that gives three constraints which are enough to impose the Wilson 
loop to be trivial, Hf = 1. Hence that proposal enforces flatness around the node, and thus 
implicitly takes into account both the scalar and vector constraints. 

Let / be a face and v a node of A on the boundary of / where ei, 62 meet. The curvature 
tensor is replaced with 

where 5/ is the Wilson loop around / based at v, and R the vector representation. The products 
of fields EfEj becomes the products of fluxes through ei and 62, and we define (see Fig. 7) 

Hvj = Ei{6ij-R{gf)^,)Ei^. 

For a face with k vertices, there are k such constraints (-ffi,„,/)n=i,...,fc around /. Since they probe 
the matrix elements of R{gf) which has three independent degrees of freedom, we expect only 
three such constraints to be independent. That was proved in [29], and the key point is that the 
differentials of these constraints around the solution = 1 is a system of rank 3. That result 
was used to prove that the constraint algebra is first class, hence it generates a gauge symmetry, 
on the 2-sphere triangulated by the boundary of a tetrahedron (the generic case being unclear 
so far). 

6.2.2 Euclidean flat geometry 

The kinematical Hilbert space gives a quantization of the intrinsic geometry. Quantizing the 
GauB law (2.1) (i.e. imposing local SU(2) invariance) is indeed quantizing the law of cosines (2.2) 
(the Euclidean formula which evaluates angles of triangles in terms of their lengths). We now 
discuss the extrinsic curvature which is actually computed from the flatness constraint. 
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Figure 8. On the left, the dashed lines represent three triangles meeting at a node, and the dual 
graph, in normal lines, depicts a face with three vertices. On the right, the geometry of the three 
triangles is pictured. Flat geometry is characterized by the fact that the dihedral angle Qtit2 between the 
triangles ti, t2 is evaluated from the three angles 4>ee' at the apex, through the Euclidean formula (6.2). 



Physically, what we expect from the classical flatness constraint F*^ = is that S with its 
geometry can be locally embedded into flat 3-space. For the triangulation dual to A, it means 
that the {3d) dihedral angles &tit2 between two adjacent triangles ti, t2 should be given as 
a standard function of the dihedral 2d angles {(pee')- A typical situation is depicted at Fig. 8, 
for a trivalent node of the triangulation. The triangles ti, t2 are dual to the vertices vi, V2 
of A, which are linked by the edge e. The classical angle between the two triangles ti, t2 can be 
computed from the three 2d angles (pee' around the node which is dual to the face /, using the 
formula 



cos et,t,{E,E) = ^""'^^1':^ COS^e,e-^^^e,e _ ^g 2) 



Sm (yeie ^'■'-^ Ye2e 



We have written it as a function of the flux variables since all angles (pee' can be evaluated from 
them, without writing down holonomies^^. 

There is another natural notion of dihedral angles which makes use of flux variables and 
holonomies. Since the fluxes represent the directions of the edges of the triangulation, the 
normal Nt to a triangle is given by the cross product of two of them. Following Fig. 8, Nl_^ = 
e^jj^XiX^_^, and N^^ = e^-^Xi^X'^. It suggests to look at the dot product of these normals as 
being the cosine of the dihedral angle between ti and t2- But for this, it is necessary to transport 
them in a common frame, say along the edge e to the vertex vi. So we define the angle 9t^t2 
as 

„ ,p ~ , Nt, • R{ge)Nt, 
cos9t,t^{E,E,g) = , • 

I ^^tl I 1^^*2 1 

This is a function of the holonomy ge since one has to compare the local embedding of ti to 
that of t2- The result can be partially evaluated in terms of the 2d dihedral angles, because it 
is a dot product of two vector products in M'^, 



coset,t2{E,E,g) 



El ■ R{ge)E2 - cos (pe^e COS (pe2e 



Obviously, one would like to this notion of dihedral angle to coincide with the standard for- 
mula (6.2). This is what the Hamiltonian does [29] 

H^j = sin (pe^e sin (pe^e (cos Gt^t2 - cos Ot^ti) ■ (6.3) 



^'^In addition, when tlie closure relation (2.1) holds, the angles are determined as functions of the lengths 
'e = \Ee\), and hence so are the angles 8*1*2 ■ 
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Figure 9. The tctrahcdral graph on the 2-sphere. 



Hence, dihedral angles are indeed computed from the 2d angles with the usual formula of flat 
Euclidean geometry. 



6.2.3 The Wheeler— DeWitt equation and the pentagon identity 

Let A be the dual graph to a triangulation of the 2-sphere by the boundary of a tetrahedron. 
It has four vertices (dual to the triangles), six links (dual to the six edges), and four triangular 
faces (dual to the nodes). A spin network function has six group elements as arguments, and 
six spins labeling the edges. Since the vertices are three- valent, the intertwiners are completely 
determined by the three spins meeting at each node, up to normalization. A common choice is 
to write the spin network function like 
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The range of summation is —je < Oei ^ je for the magnetic numbers on each link. The 
quantities into brackets are the components of the intertwiners in the usual spin and magnetic 
number basis, and are known as Wigner 3jm-symbols [112], while (j, a\g\j, b) denotes the matrix 
elements of the Wigner D-matrix with spin j. Any state ijj can be seen as a wave function and 
decomposed in that basis. 



1) 



n(2je+i) 



.e=l 



V'(jl,---,j6)Stet^(fi'l,---,ff6)- 



We now consider the node v where the edges 2 and 6 meet in the tetrahedral graph of Fig. 9 
and the face / with 1, 2, 6 on its boundary. The associated H2Q = E2 ■ Eq — E2 ■ R{gf)EQ has 
been quantized in [29]. E2 • Eq is diagonal on spin network functions, but E2 • R{gf)EQ produces 
some shifts on the spin ji due to the presence of matrix elements of gi in R{gf). The latter acts 
by multiplication on spin networks, so that shifts are due to re-coupling between the spin ji 
and the spin 1 carried by the vector representation R. In the end the Wheeler-DeWitt equation 
Hi2\ip) = is in the spin network representation 



^+i(ji)V'(ji + 1, • ■ ■ , Je) + ^o(ji)V'(ji, ■■■,36)+ ^-i(ji)V'(ji - 1, • • • , Je) 
The coefficients are A_|_(ji) = jiE(ji + 1) and A_(ji) = (ji — l)E{ji), with: 



0. 



(6.4) 



^^(Ji) = [((J2 + 33 + I? - jf) 0! - {32 - J3)') ((J5 + Je + 1)' - jl) U! - (J5 - Je)') 
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and 

Ao{j) = (2ji + I){2[i2(i2 + 1)J5(J5 + 1) + J6(J6 + l)i303 + 1) " jlUl + I)j4(j4 + 1)] 

- [J2O2 + 1) + jsUs + 1) - iiOi + 1)] [i5(j5 + 1) + jeUe + 1) - Ji(ji + !)]}• 

Hence the Wheeler-DeWitt equation is a second order difference equation, which can be solved 
from a single initial condition^®. 

This equation is sometimes taken as the definition of the Wigner 6j-symbol, ■ ■ ■ , je) = 
{ilisie}' ^ basic re-coupling coefficient of quantum angular momenta theory. If one 

wants to find a basis of the tensor product ji (S" j2 <Xi ji , one can choose a spin which 
lies in both tensor products ji (8) j2 and J4 <^ js . Another basis is obtained by choosing a spin jg 
shared by ji (g) and j2 (8) J4. The 6j-symbol gives the corresponding change of basis. 

The difference equation is thus the standard recursion relation on the 6j-symbol. The latter 
comes from the Biedenharn-EUiott (pentagon) identity [112], which is the key to re-coupling 
SU(2) angular momenta. The above result is actually the achievement of some efforts to prove 
a conjecture relating the Biedenharn-Elliott identity to the Wheeler-DeWitt equation (see [19]). 
The main idea behind the conjecture was that there exists a state-sum model, the Ponzano- 
Regge model [96], which is known to give the partition function of 2+1 gravity on a 3- manifold 
[59, 88, 90]. It is built on labeling edges of a triangulation with spins, assigning a 6j-symbol to 
each tetrahedron and summing over all labelings. When finite, it enjoys an invariance under 
some changes of the triangulation. Those changes are generated by 3-2 Pachner moves. The 
elementary move is the following. Take three tetrahedra glued along a common edge, each pair 
sharing a triangle, and remove those three triangles together with the common edge, and instead 
glue a new triangle dual to the removed edge. One ends up with only two tetrahedra glued along 
a triangle. This move preserves the topology, and the corresponding invariance of the Ponzano- 
Regge model is just due to the Biedenharn-Elliott identity. The latter thus appears as the 
algebraic translation of the Pachner move. As the 3-2 Pachner move is part of the topological 
invariance^^, a similar relation was expected in the Hamiltonian formulation. 

The relationship to flat connections on the tetrahedral graph is most conveniently seen by 
going back to wave functions of the holonomies. 



n(2je+i) 



J4 J5 J6 



31 32 33 [ lie}/ „ N 

Stet (91' • • • iffej 



It is an interesting feature that the solution is determined by one initial condition and instead of two. Looking 
at the lowest possible value of the spin ji, that is: = max(|j2 — jal, Ijs — ja\), it turns out that the lowering 
coefficient yl_i(j™'") evaluated on this spin is zero. Thus, the recurrence can be implemented starting from the 
initial value on this lowest spin. This property is important with respect to the asymptotic limit. There, it is 
known [97, 103, 112] that the 6j-symbol oscillates like the cosine of the Regge action Sn{je) ~ X]e(i<= + for 
a tetrahedron with lengths (je + |), 

fjl J2 jsl ^ 1 1 /p»(SR(j„) + f) _)_p-»(gR(je) + f)\ 

Ij4 is je / ~ ^12TvV{je) 2 V J' 

Having a second order equation, it may have been expected to be solved in the asymptotics by both the positive 
and negative exponentials of the action. But the model is fully independent from the orientation: there is only 
one required initial condition, because of the vanishing of A_i(j™'"), leading in the large spin regime to the cosine 
of the action. 

^^The full topological invariance needs the 4-1 Pachner move, but unfortunately it is divergent, because the 
model as itself does not gauge fix the non-compact gauge symmetry at the root of the topological invariance. 
Gauge-fixing in the Ponzano-Regge model is described in [59] and f initeness in [22] . 
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The above state is exactly the projector onto flat connections for the chosen graph and topolo- 

As a remark, notice that the naive state would have a delta function on each of the four faces 
of the graph. However, that would be ill-defined since one delta function would be trivially 
satisfied thanks to the other. So it is usually removed by hand. Here, the state is well-defined 
as a distribution and that is an outcome of our approach. 



6.2.4 Extra remarks 

Symmetry in the Boulatov model for 3d gravity. The latter is a group field theory, i.e. 
a non-local field theory on a group, whose Feynman expansion generates a sum over simplicial 
complexes supporting spin foam amplitudes, see [18], and more remarks on Section 6.3.3. In- 
terestingly, a global symmetry of the model has been found in [18] and there argued to be the 
translation of the diffeomorphism symmetry of 3d gravity in group field theory. The symmetry 
involves an action of the Drinfeld double of SU(2), and there is indeed a quite well established 
approach to quantize 3d gravity by going from the classical group ISU(2) (which takes into ac- 
count both the standard gauge symmetry and the translation of the co-triad) , to the quantum 
group P SU(2) at the quantum level. We do not comment further as we are not expert, but 
refer to [18] instead. 

A direct connection with the quantization of the scalar constraint we have just presented 
is obtained from the way the global symmetry of [18] is expressed in the spin basis. It is 
a recursion relation on the 6j-symbol, which can be proved from (6.4). It also derives directly 
from an operator acting on wave-functions, with the corresponding quantum constraint, 

[tr,(g/)-(2j + l)]V' = 0, 

where gj is the Wilson loop around a plaquette. In the spin basis, for a plaquette with n links 
labeled e = 1, . . . , n, it gives 
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,e=l 
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(2j + l)^{ke). (6.5) 



That equation is obviously solved by the evaluation of the spin network on the identity. 

There is an appealing geometric interpretation of that equation, as generating a quantum tent 
move evolution. We use the Ponzano-Regge interpretation [59] of the 6j-symbol: it is a weight 
associated a tetrahedron whose edge lengths are given by (je+1/2). A Ponzano-Regge amplitude 
is then obtained by taking products of 6j-symbols corresponding to a gluing of tetrahedra, and 
summing over the spins of the internal dual edges, keeping those on the boundary fixed. 

That interpretation shows that the equation (6.5) is actually computing the Ponzano-Regge 
amplitude for a piece of a 3d triangulation built as follows. The plaquette is dual to a vertex s 
of a 2d triangulation. A new triangulation, identical to the first, is obtained by evolving s into 
another vertex s* . The piece of 3d triangulation is contained between the two triangulations, 
and is made from n tetrahedra all sharing the edge {ss*). 

This evolution process in background independent approaches to lattice gravity is presented 
in [13] at the classical level, and further details are given in Section 6.3.1. In the equation (6.5), 
the evolved state ip{ke) is obtained by summing over all admissible values of the spins (je) on 
the initial surface, and the constraint just imposes that the physical state must be invariant 
under tent moves weighted by their Ponzano-Regge amplitudes. 



■^°If one gauge-fixes, using the SU(2) action at each node, the group elements gi, g2, gs to the unit, the delta 
functions impose g4 = gs — ge = 1, which is the unique flat connection up to gauge on the sphere 
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Notice that since we are deahng with a flat spacetime, the quantum length of {ss*), i.e. 
j + 1/2, is a gauge degree of freedom, which simply corresponds to the lapse. It can be chosen 
arbitrarily, and naively summing over all values of this length would result in a divergence^^. 

Thiemann's Hamiltonian in 3d gravity. The regularization and quantization of the 
scalar constraint has been performed by Thiemann [109] following the prescriptions he used in 
3+1 gravity. The construction is mathematically well-defined, as the continuum limit exists 
and there are no anomalies. The aim was, like in the new proposal we have presented, to use 
2+1 gravity as a testing ground, to be able to compare with the already known quantization. 
It turns out that the solutions to Thiemann's quantum constraints contain the usual solutions 
of [115], but also other, unusual solutions. 

Non-zero cosmological constant. The major remaining challenge for the loop quanti- 
zation of pure three-dimensional gravity is to include the cosmological constant A 7^ 0. The 
partition function and expectation values of observables are known in the case A > using the 
Chern-Simons formulation [113] which makes it clear it is still topological. In particular the 
partition function is given by the Turaev-Viro model, which looks similar to the Ponzano-Regge 
model, expect that all objects from SU(2) representation theory are changed to their counterpart 
for the quantum group C/q(SU(2)), with q = e^^. Hence, it seems that from the loop quanti- 
zation view, one has to trade off Clebsch-Gordan coefficients and 6j-symbols (and maybe spin 
networks themselves) for their g-deformed versions. It is still unclear why and how that should 
be taken care of precisely by the constraints in the loop quantization. Some recent progress have 
been obtained in [89], where the Kauffman bracket is derived by combining standard SU(2) spin 
networks with the Chern-Simons (non-commutative) connection. Interestingly, this comes out 
just from kinematical considerations. While this strengthens the robustness of the loop quanti- 
zation itself, it also leads to some puzzle. To extend the Kauffman bracket, which is defined for 
the spin 1/2, one has to tensor several spins 1/2. However, it must be followed by some sym- 
metrization process, and the usual choice in LQG is the standard symmetrizer of SU(2). But 
to get the correct extension of the Kauffman bracket, one would have to use the g-symmetrizer 
instead. There is no reason at the kinematical level to do so, and it is further unclear why the 
dynamics would change the symmetrizer. 

In the continuum, turning on A implies a deformation of the constraint algebra where the 
GauB law is unchanged but the curvature is obviously non-zero anymore, 

Kb + ^e'jj.eie'^ = 0, 

where is the pullback of the triad to the canonical surface. The algebra is first class, it is a Lie 
algebra corresponding to S0(4), where the Gaufi law generates rotations and the constraint on 
the curvature generates the Euclidean boosts. 

The most natural regularizations of the constraints in the Hilbert space of the loop quantiza- 
tion produce some anomalies in the algebra [89, 95]. It seems to indicate that since the algebra 
itself is changed, the usual regularizations (used for A = 0) must be amended to reproduce a first 
class algebra. 

A hint may come from representation theory again. When A = 0, we have two types of 
constraints, the Gaui3 law and the flatness, whose geometric contents are respectively the law 
of cosines (the formula which evaluates angles of triangles from lengths) (2.2) and a formula 
to evaluate dihedral angles between triangles from the angles within triangles (6.3). At the 
quantum level, those constraints generate difference equations solved by classical SU(2) re- 
coupling objects. It is quite clear for the closure relation: if ei, 62, 63 meet at a node, the 
constraint {Ei -\- E2 + Es)\ip) = produces recursion relations satisfied by Clebsch-Gordan 



^^For a triangular plaquette, this is nothing but the well-known divergence of the 1-4 Pachner move already 
observed by Ponzano and Regge. The quantum tent move presented here is a gauge-fixed version. 
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coefficients (or 3mj-symbols). Moreover, we have described above a regularization of the scalar 
constraint which leads to a difference equation on the 6j-symbol. 

Geometrically, the GauB law regulated as the closure relation for the triangle (2.1) only makes 
sense for a flat triangle. If geometry is homogeneously positively curved, we would expect the 
regulated GauB law to describe spherical triangles instead. In terms of lengths and angles that 
should be the spherical law of cosines, 



cos £e-) — COS le-, COS ^go 
COS0eie2 = —n . 

sm lei sm 



where ig is the spherical length of the edge e. From that view, the second constraint, that 
relating dihedral angles to angles of triangles, is actually unchanged, (6.2), 

cos ((>e^e2 - COS (/leje COS (/"eje 

cosGtita = — — . 

sm 0eieSm (pe^e 

To make those ideas more precise, in particular to quantize, it is necessary to know the expression 
of the above geometric relations in terms of holonomies and fluxes, or to have the symplectic 
structure which involves spherical lengths. 

Conjecture. Based on the results we have reviewed in the case A = 0, one can conjecture the 
following. Recursions satisfied by q-deformed Clebsch-Gordan coefficients and the q-deformed 
6 j -symbol are quantizations of the above formulae for spherical triangles and tetrahedra. 

Finding a regulated Gaufi law whose geometric content is curved triangles is certainly the 
way to get an anomaly-free regularization. This is actually what comes from the covariant 
analysis of discretized gravity by Dittrich and collaborators which we describe in Section 6.3. In 
a nutshell, the discrete action with A 7^ has no gauge symmetry when flat simplices are used. 
To restore the exact symmetry at the discrete level while turning on A > 0, it is possible to use 
a discretization based on homogeneously curved simplices. Geometrically, that is very natural 
and quite obvious: that means that translating the nodes of the lattice is a symmetry if it is 
done by moving on the curved manifold. 

Four-dimensional extension. A direct extension of the new Hamiltonian is presented 
in [27]. Its form is the same, built on holonomies and fluxes. It still acts on cycles independently, 
and thus necessarily related to the topological model which imposes flatness of the connection. 
On the boundary of a 4-simplex, triangulating the 3-sphere, the Wheeler-DeWitt equation is 
a difference equation solved by Wigner 15j-symbols. 

More generally it is easily seen that the Hamiltonian we have discussed produces difference 
equations which are solved by Wigner symbols. It is actually not surprising, but it is instead 
a confirmation that the formalism is under control. Indeed, Wigner symbols are evaluation of 
spin networks on the trivial connection, and as such they naturally arise in the quantization of 
the BF model. 

Large spin asymptotics and loop quantum gravity geometries. The large spin limit 
of the recursion (6.4) makes the geometric picture quite clear, 

[Ae + 2(l-COSee)]V(il,-.-,j6) =0, (6.6) 

where Ag is the discrete Laplacian on the edge e, and Qe the dihedral angle at e, computed 
from the lengths. Such an equation can be straightforwardly derived from a simple dynamics 
in Regge calculus. The canonical variables are the edge lengths Ig and momenta Og, and the 
constraint reads 

COS^e — COS0e(/) = 0, 
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which means that the momenta are fully determined as dihedral angles between triangles for 
flat space. 

A similar constraint was proposed for flat space in 3+1 dimensions in [47], with lengths 
replaced with triangle areas, and it was argued that they form a first class system. As a con- 
vergence of ideas and approaches, it turns out that the BF model provides a realization at 
the quantum level of such a constraint for flat space. Indeed, the large spin behavior of the 
difference equation coming from H^j on the boundary of the 4-simplex, using coherent state 
methods, takes exactly the form (6.6), where spins label triangles and geometrically stand for 
quantum areas instead of quantum lengths [27] . 

However, we would like to stress that the Hamiltonian on holonomies and fluxes goes beyond 
Regge calculus in the 3+1 dimensional case, in a sense we now explain. In three dimensions, 
all flux configurations which satisfy the closure equation are Regge geometries, determined by 
length variables. We simply say that loop quantum gravity geometries are Regge. But the situ- 
ation is different in four dimensions, because a tetrahedron is not determined by its areas. The 
result of the full analysis [47] is that SU(2) loop quantum gravity geometries are discontinuous: 
each tetrahedron unambiguously describes a flat Regge geometry, but different tetrahedra will 
generically describe their common triangles with different lengths. Such geometries were later 
called twisted geometries [61]. 

Instead of going further on the kinematical aspects (see [47, 100]), we simply want to empha- 
size that, as far as flat space is concerned, the Hamiltonian is well-defined on the whole phase 
space of loop quantum gravity, while the large spin behavior (6.6) only describes the geometric 
sector, where geometries are really of the Regge type. For non-Regge configurations, the exact 
difference equation takes a different asymptotic form, which is consistent with asymptotics of 
different symbols in that regime [20]. In the geometric sector, (6.6) can be taken as a criterion 
for a quantum geometry model to be approximated with some quantum Regge calculus in the 
large spin limit. 

More on recursion relations. The study of H^j and other recursion relations on Wigner 
symbols has suggested a nice feature of quantum BF theory, that geometric properties of quan- 
tum flat space are encoded into recursion relations coming from group representation theory. An 
example, discussed in [27, 31, 33] is the closure of the simplex. In the asymptotics, the idea is to 
combine (6.6) with the classical closure which reads det(cos 0e(/)) = 0, i.e. the vanishing of the 
Gram matrix. Remarkably, the lOj-symbol at the core of the four-dimensional Barrett-Crane 
model is known to satisfy such an equation, but in an exact way [33]. This leads to expect that 
in the geometric sector of full quantum gravity, geometric properties are still probed through 
difference equations which arise from representation theory. 

6.2.5 What we learn 

In addition to its own interest as a new way of quantizing 2+1 gravity using the scalar Hamil- 
tonian, this analysis gives the following insights. 

• The spin 1 is the natural spin to regularize the curvature. It makes it possible to solve the 
Wheeler-DeWitt equation for all values of the state on integers. Hence, that solves the 
model for SO (3). By contrast, working with a spin 2 regularization produces a difference 
equation with shifts ±1, ±2, so that more initial conditions are needed to get the state on 
all integers [29]. Solving the model for SU(2) instead of S0(3) however requires a spin 1/2 
regularization of the curvature, which is accessible in the spinor formalism [84]. 

• The splitting between spatial diffeomorphisms and the scalar component finds an unex- 
pected realization in this exact model. Indeed, the usual tension is that the curvature at 
a node of a triangulation is regularized by Wilson loop around it, while the scalar Hamil- 
tonian is defined by contracting the curvature with the normal to the surface, H ^ n ■ F . 
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The way-out which turns out to be successful is to take a different normal for each triangle 
meeting at a node, three of them being linearly independent. From that view the above 
construction is more something which mimics the scalar constraint on spin networks. 

• The asymptotic equation (6.6) is a key feature of the quantum simplex in the large quantum 
number regime, and can used as a criterion for a model to reproduce Regge calculus there. 

• Classical geometric properties of flat space can be implemented at the quantum level 
as difference equations coming from group representation theory. This strengthens the 
core of spin foam models where spacetime dynamics is expressed with objects from group 
representations. 

• The loop quantization supports a non-trivial Hamiltonian, not only in the asymptotics but 
also for twisted geometries in four dimensions. 

• No divergences are present, because the generator of the gauge symmetry is quantized and 
not the projector. That improves the spherical case with respect to [88], and it should be 
a good framework to study perturbations with respect to the topological order. 



6.3 Beyond the topological model 
6.3.1 Broken symmetries 

The above BF model is topological and enjoys an exact lattice realization which allows the 
evaluation of transition amplitudes between spin network states using a single cell decomposition 
of spacetime (containing the graphs). But generically in a theory which does have local degrees 
of freedom, a lattice formulation is a truncation which misses some degrees of freedom and hence 
breaks the symmetries. In a series a papers Bahr and Dittrich and collaborators [13, 14, 15, 17], 
and Dittrich and Hohn [46], studied the discretization of reparametrization invariant systems, 
especially the Regge approximation to general relativity. Main issues are how to restore broken 
symmetries, and how to define a canonical discrete formulation which reproduces the dynamics 
of the discrete action. Here we sum up some ideas of [13]. 

Starting from a triangulation of the initial canonical hypersurface, a local evolution with 
discrete time steps can be defined in terms of tent moves. Let A(^^_i^ be the initial triangulation, 
and s a node. A tent move on s creates a new node s* , and new links which join s* to all the 
nodes s was connected to. That defines a new canonical triangulation A^^_-|^^, on which the 
number of variables is unchanged. Equivalently, the move glues a piece T of a d-dimensional 
triangulation between A(^_i) and A*^_-^y and the link {ss*) is called the tent pole. 

The independent variables on T are the lengths. To get momenta and equations of motion 
on A(^^_i-j and A^^_-^^, one uses the action on T as a generating function. When d = 3, it reads 

St = - - ^-'i'^e + 

eCintT e(ZdT v 

The deficit angles around the edges e in the bulk of T and on the boundary are 
ee= f 27r - ^ i?ei, j , and = f tt - ^ i^e-u j , 

where "Q^v is the dihedral angle hinged at e in the tetrahedron v. Vy\s, the volume of and A 
the cosmological constant. Denote Sj the nodes connected to s, with lengths li and /* the tent 
pole length. The momenta are then 

__5SV ^_dST OSt 

^'~~~dk' ^'~~dl[' ^~~~dF' 
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From the equations of motion, we get p = 0. Moreover, performing a second tent move, it can 
be seen that there is a single well-defined momentum per edge. 

Symmetries on a solution can be observed from the existence of zero eigenvalues of the Hessian 
matrix evaluated on that solution. In the case of vanishing cosmological constant, a symmetry is 
found in [13] around the flat solution. It corresponds to moving a internal vertex in flat 3-space 
(after several tent moves are performed), and those degrees of freedom are interpreted as the 
lapse and the shift. That translates into constraints on the canonical data on A(^_x) only. 

Pi = --n- + T^ss,{h,hj), (6.7) 

in agreement with the fact that 3d gravity supports a discretization with exact symmetries. 
Those constraints state that the momenta are identified with the dihedral angles determined by 
the lengths. This is analogous to the constraints (6.3) considered in the previous section. We 
thus see a convergence of idea coming from different approaches to the same issue. 

However, when turning on the cosmological constant, it was observed that there are no 
symmetry around homogeneously curved solutions. Therefore, there are no constraints anymore, 
but instead pseudo-constraints. That means that the equations of motion from T translates in 
the canonical setting into equations which depend on the data living on two consecutive time 
slices, and not a single one like in (6.7). 

6.3.2 Restoring symmetries and the continuum limit 

In that example (homogeneous curvature), it is nevertheless quite clear how one can discretize 
in an exact way. The correct symmetry should be translation of a vertex not in flat space but 
following the curved solution. Hence, one should use curved simplices instead of the flat ones, 
so as to get the constraint 

where the dihedral angle "di^} is evaluated within the curved tetrahedron. 

The main result of this approach is a canonical framework which reproduces the discrete 
action. The presence of symmetries is therefore determined by the choice of the discrete action. 
That view was further expanded in [46] by studying the expansion of Regge calculus around 
a background metric, and in [14]. In the latter, the key idea is to restore the symmetries which 
are broken upon discretization by going to the perfect action. That perfect action is reached, at 
least formally, in the continuum limit. Hence, it can be approached by a coarse-graining process 
at the classical level. Let A2 be a triangulation whose edges are denoted E, and Ai a finer 
triangulation with edges e. We will obtained from the Regge action on Ai with flat simplices, 
an improved action on A2. Let us fix the lengths Le, and solve the equations of motion for Ig 
with a homogeneously curved geometry determined by A. The action on the coarse-grained 
triangulation will be obtained thanks to the constraint Yl = Le- The improved action is 

eCE 

then 

eCE 

where the lengths /* are solutions of the curved geometry. The action for A2 consequently 
takes into account the dynamics of the finer triangulation. Repeating the process many times, 
the internal geometry of the simplices of the most coarse triangulation receives more and more 
contributions from the curvature of the finer ones, and thus in the appropriate limit, each 
simplex becomes homogeneously curved. That leads to the Regge discretization which has the 
exact symmetries. 
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While the topological BF model supports an exact discrete version, gauge symmetries are 
generically broken. Restoring them (and hence getting the correct number of degrees of freedom, 
in particular of propagating degrees of freedom) necessitates to go to the continuous limit. It 
can be reached as above by coarse-graining (recursive integration of the equations of motion) 
to yield at least formally a discrete action with exact symmetries, or by taking the number of 
simplices to infinity. 

The work of Dittrich and collaborators was recently re-interpreted by Rovelli [99]. The 
desired reparametrization invariance (diffeomorphisms) to be recovered was there called Ditt- 
invariance. One key feature is that the continuum limit is reached by sending the number of 
discrete steps of the parameters of the theory to infinity, while the physical time step is sent to 
zero in a dynamical way. In the case of the harmonic oscillator, portions where the trajectory 
can be approximated by a free particle need very little time steps. On the contrary portions 
where the potential is important need more points to describe the trajectory accurately. It turns 
out that the discretization in the arbitrary time parameter together with energy conservation 
leads to that picture in terms of the physical time parameter. 

The discretization and coarse-graining process is discussed for reparametrization invariant 
systems at the quantum level in [17]. 

6.3.3 The continuum limit through matrix and tensor models 

The notion of continuum limit at the quantum (or statistical) level has been given a sense in 
matrix models (for two-dimensional surfaces) [44] and tensor models (in higher dimensions) 
[4, 67, 102]. These are field theories which generate sums over triangulations, and as such the 
latter can be controlled. The free energy of a typical model reads 

F = T. E s{Tn)A{Tn). 

n>0 triangulations T„ 

with n simplices 

Here T„ are triangulations of the sphere^^ with n simplices, s(T„) is a combinatorial factor, 
and A{Tn) the amplitude assigned by the model to a given triangulation. The feature we 
want to stress here comes from the chemical potential When it is very large, only small n 
contribute. However, when it decreases, the sum loses its summability. This is the regime 
where F is dominated by an infinite number of simplices, and it is obtained for a (non-zero) 
critical value of ^ — )• /Xc. This continuum limit is approached like F ~ (e"^'' — e~^) . In the 
simplest models, where A{Tn) = 1, the sum at fixed n just counts the number of triangulations 
with a given number of vertices, so that 7 is known as the entropy exponent. In the matrix 
model, it is 7 = —1/2 and in higher dimensions, 7 = 1/2. 

Such models have a natural interpretation as dynamical triangulations, in the regime of small 
Newton constant and large cosmological constant. It is also possible to modify the ampli- 
tudes A{G) to accommodate loop quantum gravity, or rather spin foam, amplitudes. That is 
the purpose of group field theories, which re-sums spin foam amplitudes, and much of the recent 
progress in tensor models have been motivated by the group field theory program. 

7 Conclusion 

Since the original definition of Hamiltonian constraint by Thiemann, several new approaches to 
tackling quantum dynamics in LQG have emerged, and the quantum constraint, still the most 

■^■^They are isolated from the full sums over triangulations of all topologies using the large A'^ limit [42, 68, 74], 
where A'^ is the typical size of the matrices or tensors. The triangulations to be summed over may moreover be 
a specific subset of all triangulations of the sphere [30]. 
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popular of all its rivals in mini-superspace models, has taken a backseat in the full theory. The 
infinite number of ambiguities, lack of spacetime covariance at quantum level, and ultra-locality 
are some prime reasons why it is even now considered to be a highly unsatisfactory definition. 

We have argued that the quantum spacetime covariance, i.e. the off-shell closure of the 
constraint algebra, could be achieved if one was ready to broaden one's perspective on all 
possible choices that are available in defining composite operators, thanks to the discontinuity 
of representations and its inherent non-local nature. Already in loop quantum cosmology and in 
parametrized field theory, one sees the fruits of such endeavours and ends up with regularization 
schemes which remarkably enough have never been considered in the full theory. 

Applying these lessons to the diffeomorphism constraint bears some real fruits. A precise 
definition of the quantum constraint which generates infinitesimal diffeomorphisms on states 
in the LM habitat, and whose kernel matches with the well-known diffeomorphism invariant 
space "Hdiff- Once again the choices one has to make in order to get such an operator are very 
different then the ones which have been looked at in the literature so far. 

We have then presented a different road based on the loop quantization of the BF topological 
model. It includes pure 2+1 gravity, where a recent Hamiltonian, [27, 29], was found to mimic 
the scalar constraint while reproducing the standard results. That shows that the framework 
supports non-trivial operators, with a nice geometric content: the result is a quantization of flat 
(Euclidean) geometry. More generally, that model shows that classical properties of flat geome- 
try take the form of difference equations from group representation theory, in the spin network 
quantization scheme. It also shed some light on the quantization ambiguities, such as curvature 
regularization and its spin, and is close to Smolin's proposal [104] to avoid ultra-locality. How- 
ever, to go further, the continuum limit should be investigated, beyond the topological model. 
We have presented some directions in the context of discrete gravity, where symmetries are 
broken, but can be restored upon coarse-graining, leading to the perfect action. 

The remaining issues in 2+1 gravity are mainly to write explicitly the quantum algebra 
of the new Hamiltonian. It should be interesting since it does not implement the splitting 
between spatial diffeomorphisms and the scalar component in the usually expected way. That 
was a surprise and it should certainly be considered seriously also in 3+1 general relativity. 

This unexpected splitting of the algebra observed in 2+1 gravity speaks in favor of the spin 
foam formalism, because it is a space-time formulation on cell decompositions, and hence natu- 
rally compatible with the geometry arising from the new Hamiltonian. The difference equations 
coming from group representation to describe quantum flat geometry also point towards spin 
foam models, since they are built from group Fourier transforms. In particular, the asymptotic 
form of the Wheeler-DeWitt equation derived through the new Hamiltonian (for the topolog- 
ical model) is satisfied on the quantum 4-simplex by some spin foam models [27]. However, 
in 3+1 dimensions, exhausting the gauge symmetries is subtler (due to reducibility). Hence, 
in spite of those quantum difference equations which match with the spin foam model for 3+1 
BF theory [91], the latter is not well-defined yet: it is not known how to gauge-fix it properly. 
We think it is a serious issue which should be resolved if one wants to understand what the 
diffeomorphism symmetry becomes in the flat or almost flat sector of the full theory, and how 
the internal gauge symmetry of BF theory is broken on non-flat solutions. 
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